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Abstract. The auto-covariance function of the differences between the observed and predicted ecliptic longitudes of 
‘upiter indicates that these residuals contain a trigonometric term of period 4520 days, compared with 4333 days for 
‘upiter’ s period, and amplitude 0’25. No currently known source of systematic error would account for such a period. 

The residuals for the other outer planets contain no period as marked as the 4520 day period in Jupiter. A secular increase 
n the eccentricity of 1/3 X 107° per century, together with a motion of the perihelion in the direction of Jupiter’s revo- 
ution about the sun and amounting to 5” per century will, coupled with .a term whose frequency is equal to Jupiter’s 
inean motion, produce such a period. If such secular perturbations were produced by a constant anomalous acceleration 
bf Jupiter, the magnitude of the acceleration required would be about 1078 times the Jovian centripetal acceleration. 

iThe data for the other outer planets are inadequate to either confirm or deny the existence of such an acceleration in 


(heir motion. 


| There appears to be considerable misunder- 
btanding on the part of physicists of the nature 
‘and degree of the observational support of gravi- 
tational theory. For example, it appears to be 
commonly believed that the observations of plan- 
etary motion agree with computed orbits to the 
‘accuracy of the observations. On the other hand, 
it has long been known by the astronomers that 
ithere are sizable systematical discrepancies be- 
‘tween computed and observed orbits (Clemence 
/1951). It has been assumed that these discrep- 
vancies have their origin in faulty perturbation 
‘calculations and systematic observational errors 
‘of unknown origin. 

It thus appears that physicists have great con- 
fidence in gravitational theory, in large part, 
because of an erroneously presumed perfect agree- 
ent between astronomical observations and 
heory. On the other hand it is likely that the 
confidence of the astronomers is, at least in part, 
based on the confidence of the physicists respon- 
sible for the basic theories. 

Till recently, the theoretical orbits with which 
the observations were compared had been ob- 
tained by lengthy perturbation calculations which 
gave the position of the planets less accurately 


than they were observed. For the five outer 
planets, these orbits have now been obtained by 
numerical integration on an electronic computer 
by Eckert, Brouwer and Clemence (1951). This 
integration procedure presumably reduces any 
uncertainty in the theoretical predictions of the 
positions of the outer planets to much less than 
the uncertainty in the observations. The calcu- 
lation was checked by computing for all times 
the total energy and angular momentum of the 
planets (Clemence and Brouwer 1955). 

The initial positions and velocities of the 
planets were first taken from existing theories, 
and the positions at convenient earlier times then 
calculated. The initial conditions were adjusted 
to improve agreement with the observations at 
that earlier time. This procedure was applied to 
Jupiter, Saturn, Uranus and Neptune, but there 
are not enough older observations of Pluto that 
the method could be applied to this farthest 
planet. The motion of Pluto given by the numeri- 
cal integration agrees with modern observations, 
but is expected to disagree with observations by 
the year 2000. The effect of Pluto on the other 
planets is so small that this should not cause an 
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appreciable error in the predicted motions of 
Jupiter, Saturn, Uranus and Neptune. 

In computing the solar attractions for the 
numerical integration, the total mass of the inner 
planets was added to that of the sun. The small 
additional effect that the inner planets have on 
the outer ones can be taken into account by a 
supplementary perturbation calculation (Clem- 
ence 1954). 

The results of the numerical integration are 
the rectangular heliocentric equatorial coordi- 
nates of the five outer planets given at 4o day 
intervals from 1653 to 2060 to g places of deci- 
mals, the unit of length being the semimajor 
axis of the earth’s orbit. 

Eckert, Brouwer and Clemence (1951), Tables 
I, 2, 3, and 4, give a comparison between obser- 
vations and the numerical integration for Jupi- 
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Figure 1. Residuals in ecliptic longitude for Jupiter and auto-covariance of residuals. : 
Figure 2. Residuals in ecliptic longitude for Saturn and auto-covariance of residuals. ai 
Figure 3. Residuals in ecliptic longitude for Uranus and auto-covariance of residuals. aH 
Figure 4. Residuals in ecliptic longitude for Neptune and auto-covariance of residuals. < 
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ter, Saturn, Uranus and Neptune. The residua 
quoted there are plotted here in the upper part 
of Figures 1, 2, 3, and 4. These residuals 
differences between observed and predicted ecliy 
tic longitudes, in the sense Observation — Ca, 
culation. . 
The residuals contain both random and sy; 
tematic components. The systematic errors al 
presumably due to systematic errors in the obse’ 
vations, to small inaccuracies in adjusting t 
initial conditions of the numerical integratioi 
and possibly to inadequacy of the theory. TE 
general scatter of the residuals for Uranus ar 
Neptune decreases with time which is to be e: 
pected as techniques improve. Such a decrea 
is not evident for Jupiter and Saturn. This su; 
gests that the residuals for the latter two plane! 
are largely systematic. 
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TABLE I, RESIDUALS FOR JUPITER 


59 June 
Residual in 
Longitude 

‘ulian Day Year Oo- Julian Day 
:7 1000.5 1779 +0o"1 239 1000.5 
,7 2000.5 1782 0.1 239 2000.5 
57 3000.5 1784 +0.4 239 3000.5 
7 4000.5 1787 OnE 239 4000.5 
$7 5000.5 1790 —0.8 239 5000.5 
7 6000.5 1793 —I.1 239 6000.5 
7 7000.5 1795 —0:1 239 7000.5 
7 8000.5 1798 0.0 239 8000.5 
7 9000.5 1801 — O35 239 9000.5 
8 0000.5 1804 —0.9 240 0000.5 
8 1000.5 1806 —0.9 240 1000.5 
'8 2000.5 1809 —0.4 240 2000.5 
\8 3000.5 1812 —0.3 240 3000.5 
\8 4000.5 1815 —0.7 240 4000.5 
\8 5000.5 1817 —0.7 240 5000.5 
‘8 6000.5 1820 —0.4 240 6000.5 
\8 7000.5 1823 —0.5 240 7000.5 
8 8000.5 1826 —0.6 240 8000.5 
8 9000.5 1828 —0.7 240 9000.5 
9 0000.5 1831 +0.2 241 0000.5 


This can be checked by looking at the signs of 
ie residuals. The residuals for Jupiter, as given 
y Eckert, Brouwer and Clemence (1951) in 
able I, are reproduced here in Table I. 

There are 60 residuals in Table I, 28 positive 
id 32 negative. If these 28 plus signs and 32 
inus signs were arranged at random, each 
rangement as probable as any other, then the 
‘obability of there being exactly k runs (k = 1, 
3, °°) can be computed and is shown in Fig- 
e 5. A run is a sequence of like signs, with a 
ro being considered to have the same sign as 
e preceding residual. The number of runs is 
us one more than the number of sign changes. 


DISTRIBUTION OF RUNS 
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jure 5. Probability distribution for number of runs in a 
random distribution of 28 + signs and 32 — signs. 


Residual in Residual in 


Longitude Longitude 
Year - Julian Day Year - 
1834 +027 241 1000.5 1888 —0'4 
1836 +0.3 241 2000.5 1891 —0.5 
1839 0.0 241 3000.5 1894 —0O.1 
1842 +0.1 241 4000.5 1897 0.0 
1845 0.0 241 5000.5 1899 0.0 
1847 —0.3 241 6000.5 1902 +0.2 
1850 —0.4 241 7000.5 1905 +0.3 
1853 —0.6 241 8000.5 1908 +0.4 
1856 0.2 241 9000.5 1910 +0.3 
1858 +0.4 242 0000.5 1913 +0.3 
1861 +0.2 242 1000.5 1916 Ox 
1864 +0.3 242 2000.5 1919 +0.1 
1867 —0.1 242 3000.5 1921 +0.2 
1869 +0.3 242 4000.5 1924 +0.3 
1872 +0.4 242 5000.5 1927 —0.7 
1875 +0.2 242 6000.5 1930 One 
1878 +o.1 242 7000.5 1932 —0.2 
1880 +o.1 242 8000.5 1935 =0.3 
1883 0.0 242 9000.5 1938 —0.5 
1886 —0.2 243 0000.5 1941 —0.5 


The most probable number of runs is 30 or 31, 
and the probability that a random arrangement 
would give less than 18 runs or more than 42 is 
negligible. In the residuals of Table I, there are 
only 10 runs. The probability of this occurrence 
is I0—8, which implies that the residuals are corre- 
lated ; they change sign less often than would be 
expected if they were random. 

One way of testing the residuals to see how 
successive residuals are correlated, and in par- 
ticular if any frequencies are present with large 
amplitude, is to compute the auto-covariance 
function s;, defined as follows: 


rt n—k 


ie kb ye VV itke 


i=1 


(1) 


St. aa 


They,,2 = 1 --- mare the residuals (with x = 60 
for Jupiter). The fourier transform of the s’s is 
essentially the square of the fourier transform of 
the 7’s, so that strong frequency components are 
enhanced relative to weak ones. If the residuals 
r; fluctuated about zero (with mean square a?) 
then, though s) would have to be positive (and 
= o”) the succeeding s;, k > 0, would fluctuate 
about zero, the root mean square of s; (for fixed 
k) being 


fig 


Sn (SW EON 12) 

If the original residuals 7; were a trigonometric 
term A cos (w? + 6) plus a random variable, then 
the auto-covariance s, would have a trigono- 
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metric term 4A? cos wi superimposed on the ran- 
dom scatter. 

The results of such an auto-covariance analysis 
are shown in the lower parts of Figures 1, 2, 3, 4. 
The residuals given by Eckert, Brouwer and 
Clemence (1951) in Tables 1, 2 (Jupiter and 
Saturn) have equal time intervals of 1000 days 
for Jupiter, 1920 days for Saturn. The residuals 
for Uranus and Neptune are not at equal time 
intervals, but are in groups of I, 2 or 3 about a 
year apart. In computing the auto-covariance of 
Uranus and Neptune, the residuals in a group 
were averaged and plotted exactly one year apart. 
A residual usually had to be moved about 0.2 
years, which should have negligible effect on any 
trigonometric terms with period more than sev- 
eral years. 

If the residuals were normally distributed with 
mean 0 and mean square o” = So, then the amount 
of scatter to be expected is indicated roughly by 
the two horizontal lines which are drawn a dis- 
tance o?/~n above and below the x axis. The 
only auto-covariance plot which shows a regular, 
definite oscillation is that for the planet Jupiter. 
There are hints of periodicities in the other auto- 
covariance plots and a fourier analysis of these 
graphs would probably exhibit them, but for 
Jupiter the oscillation is well defined and sus- 
tained for many periods. This oscillation corre- 
sponds to a term 


” F 360° ee ° 
0725 sin =, 7627) (2) 


in the residuals themselves, the phase having 
been determined by the method of least squares. 
tis the time in days from Julian Day 2371000.5 
(the time of the first entry in Table I). 

The reality of the term (2) is supported by 
Figure 6, in which is compared the auto-covari- 
ance of the residuals from Table | (the solid 
curve is drawn through the s; so computed) with 
the auto-covariance (solid circles) computed from 
these residuals after they had been permuted 
randomly. The solid circles show no regular oscil- 
lation, and appear to scatter irregularly. The 
auto-covariance of the Jovian residuals after the 
term (2) had been subtracted from them is shown 
in Figure 7. The regular oscillation is gone, and 
a few periods of another, slower, oscillation are 
left. There is still little random scatter, which 


suggests that the residuals after subtraction of | 


(2) are still largely systematic. Subtraction of 
the term (2) decreased the mean square residual 
from 0.175 to 0.148, a reduction of 15 per cent. 
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. Figure 8. Auto-covariance of Jupiter’s residuals before an 


The oscillation (2) present in the residuals o 


Table I has a period (12.4 years) somewha 


greater than the period of Jupiter itself (114 
years) in its motion around the sun. The diffe 
ence between these two periods would appear by 


be real, because there is no trigonometric tern 
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Figure 6. To show the effect that randomizing the order, 
Jupiter’s residuals has on the auto-covariance. ; 
Figure 7. Auto-covariance of Jupiter’s con before an 


after subtraction of the term 0725 sin CS — 76°7 }. 


4520 / 


: : 60° * 
after subtracting a term —0’13 sin ( SPO fates “a 
| 3 4333 + 233°} 


59 June 


ith period 11.9 years which can be subtracted 
om the residuals of Table I to remove the 
cillation present in the auto-covariance of Fig- 
re 1, The best, in the sense of minimizing the 
ean square residual, that can be done with a 
igonometric term of period 11.9 years is shown 

Figure 8. Here the solid curve is again the 
ito-covariance calculated from the residuals of 
able I, while the dotted curve is drawn through 
ie s, calculated from these residuals after a 


Tm 
° 


6 
—o'13 sin & . t+ 233°) (3) 


4333 
1d been subtracted from them (the period of 
ipiter is 4333 days). ¢ in (3) is reckoned in days 
om Julian Day 2371000.5, as was the ¢ in (2). 
he amplitude and phase of the term (3) were 
1osen to minimize the mean square residual. 
ibtraction of this term reduces the mean square 
ily from .175 to .166, and fails to remove the 
cillation in the solid curve. 

Discussion. Assuming that the residuals in 
able I have the statistically significant term 
») in them, what does the presence of such a 
rm mean? 

It is difficult to think of a source of systematic 
ror which would lead to an inequality in the 
siduals having a period just a bit longer than 
ipiter’s period. Errors in the positions of the 
andard stars with respect to which Jupiter is 
easured would lead to either Jupiter’s fre- 
iency or the difference between the earth’s 
bital frequency and Jupiter’s frequency, as 
ould systematic errors dependent upon the time 
year on the earth. Errors in the initial condi- 
ons used in the numerical integration would 
so lead to a term in the residuals having Jupi- 
r’s frequency. The inner planets can affect the 
ngitude of Jupiter by several tenths of a second 
arc, but the frequencies of the terms so intro- 
iced into Jupiter’s longitude are all much too 
gh. 

There are some theoretical corrections which 
ould be applied to Jupiter’s calculated orbit. 
ne of these is due to the fact that the force 
sting on the center of gravity of Jupiter and its 
itellites is not exactly the derivative of a gravi- 
tional potential evaluated at the center of 
avity, but contains some correction terms be- 
use Jupiter and its satellites are an extended 
stem. Another correction is for the fact that 
le perturbing planets (principally Saturn) are 
liptical in shape, not spherical. Thirdly, the 
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observed center of Jupiter is displaced from the 
center of gravity of Jupiter and its satellites, and 
it is the motion of this center of gravity that is 
given by the numerical integration. All these 
effects are negligible (Hertz 1953), and do not 
explain the observed inequality (2). 

The frequency of the observed inequality (2) 
is 286’ per day. The only frequency present in 
the solar system which could combine with Jupi- 
ter’s mean motion, 299” per day, to give a beat 
frequency of 286” per day is Pluto’s mean mo- 
tion, 14’ per day. The perturbation caused by 
Pluto in Jupiter’s ecliptic longitude is approxi- 
mately given by the formula 


07007 sin (J’ — 1) — o“002sin2(l’ —1)+--- 


(4) 


where /’ is Pluto’s longitude and / is Jupiter’s 
longitude. (The amplitudes were calculated from 
the formula given by Tisserand (1889, Vol. 1, 
p. 365), neglecting terms proportional to the first 
and higher powers of the eccentricity, and the 
inclination of the orbits.) A planet with a mass 
35 times greater than Pluto’s mass and at the 
same distance from the sun as Pluto would pro- 
duce the term (2) in Jupiter’s longitude, but 
would also produce terms of amplitude 1” in 
Saturn’s longitude and 16” in Uranus’ longitude. 
The absence of such a term in the residuals of 
Uranus (Figure 3) rules out such a planet. 

The observations of Jupiter on which the re- 
siduals in Table I are based were not made 
continuously throughout the year, but rather in 
groups spaced about 1.092 years apart (i.e., when 
the sun, earth and Jupiter were in line). Thus 
there is a possibility that the observed oscillation 
in Jupiter’s residuals is a distributed effect, and 
really represents an error of much higher fre- 
quency, the observed frequency of 1/12.4 yr 
being a beat between the ‘‘true error frequency” 
and the frequency, 1/1.09 yr—!, at which obser- 
vations are made. This “‘true error frequency” 
would have to be less than the mean motion of 
the earth by about 14” per day, which corre- 
sponds to a period of 366.5 days. The authors 
know of no factor which would lead to a periodic 
error in Jupiter’s longitude with a period of 
366.5 days. 

To get a feeling for what the discrepancy (2) 
means for Jupiter’s orbit, we can ask what kind 
of deformation in Jupiter’s orbit would lead to 
such a term. Consider a planet moving on an 
ellipse whose eccentricity is e and whose peri- 
helion makes an angle w with some fixed direction 
in space; when the planet’s mean longitude is /, 
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its true longitude v will be given by the relation 
=1+2esin(@l—w)+::: (5) 


if the eccentricity is small and the longitudes are 
measured with respect to the same fixed direction 
in space. An error in e and @ will lead to an error 


dv = 26e sin (| — w) — 2e6w cos (J — w) 
+ higher terms (6) 


in the true longitude. If 6e and 6m are constants, 
then 6y will vary harmonically with frequency n, 
where n is the mean motion of the planet, but if 
de and eém are also time dependent, for example 
de = cos wt and edw = sin wf, then the error will 
have a frequency 1 — w. This is what happens 
in the perturbation of Jupiter by Pluto; the 
angular frequency w is about 13” per day. 

This difference frequency w is small enough 
that over the 160 years covered by Table I, 
sin wt can be approximated reasonably well by 
wt, and cos wt by 1, if t is measured from the 
midpoint of the interval. Although taking 6e and 
eém in (6) to be linear functions of the time, and 
using this approximation, will not lead exactly 
to the term (2), this term can be represented 
quite well if a suitable trigonometric function 
varying with frequency 7 is added to (6). The 
secular perturbations required are an increase of 
1/3 X 10-8 parts per century in e and of 1.3 X 

’ Io-6 radians per century in edw. The latter figure 
corresponds to a positive motion of the peri- 
helion amounting to 5”’ per century. If the'secu- 
lar perturbations are taken to vanish at the mid- 
point of the interval, the term with frequency n 
which must be added to (6) is 0725 sin (w + 15°), 
where w is Jupiter’s position angle measured 
relative to perihelion. 

The validity of replacing the sine and cosine 
by the first terms in their Taylor expansions over 
the time interval in question is demonstrated in 
Figure 9, which shows the auto-covariance of a 
fictitious set of residuals which have -been con- 
structed by adding the term 


.065 ¢ sin (1.50 ¢) + cos (1.50 2) (7) 


to a set of random numbers drawn from a gaus- 
sian population. In (7), £ is measured in units of 
1000 days from the midpoint of the 160 year 
interval covered by the residuals of Table I, 
i.e., from Julian Day 2401000.5. The numbers 
1.50 and 0.065 are Jupiter’s mean motion and the 
difference frequency 13/’/day expressed in ra- 
dians/1000 days. The random numbers were 
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64, No. 127 


e dotted Iine* random numbers ' 
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Figure 9. Upper part: ‘Residuals’ constructed by add 
ing .065 sin 1.50% + cos1.50% (dotted line) to a set c 
random numbers. Bottom part: auto-covariance of thes 
“‘residuals.’’ The frequency of the oscillation is 1.44 radians 
1000 days which is to be compared with 1.50 — 0.65 = 
1.435 radians/1000 days. 


normalized so as to have the mean square fiv 
times greater than the mean square of the tern 
(7); the mean square of the term (2) in Jupiter’ 
residuals is about 1/5 of the mean square g 
sidual. Although 0.065 ¢ varies from —2 to +% 
radians, as far as the auto-covariance is con 
cerned the approximation in which (7) is r 
placed by a single cosine is not a bad one, for th 
oscillation in Figure 9 appears well defined ant 
has a frequency (1.50—0.065) radians/1000 days 

If such secular perturbations really exist (due 
possibly to incorrect planetary mass assignments 
an error in the general precession, or conceivabh 
even an anomalous, non-Newtonian, force) th 
source of the “suitable term of frequency #7 
might lie in the adjustment of the initial cond 
tions for the numerical integration. A term Like 
0.065 / sin (1.504) arising from secular change: 
in the orbit would lead to residuals whose ampli 
tude varied with time; addition of a term like 
cos (1.50 t) (by adjusting e and w) would ma 
the residuals more uniform in time withouw! 
changing appreciably their mean square. Whethel 
this is a reasonable explanation or not depends 
on the details of how the initial conditions for 
the numerical integration were chosen. 

To get an idea of how large a force is requiret 
to produce a term like (2) in the residuals @ 


produce a secular motion of the perihelion such 
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; not statistically significant). The rates of 
hange of the eccentricity and position of the 
erihelion are related to perturbing accelerations 
J and TJ normal and parallel to the velocity v 
aspectively through the formulae, valid forsmall in adjusting the initial conditions. 
ccentricity, Conclusion. The auto-covariance function of 
eo oN oT the differences between the observed and pre- 
mag Ceo) = > 60s (i—w) (8) dicted ecliptic longitudes of Jupiter indicates that 
these residuals contain a trigonometric term of 
nd period 4520 days, compared with 4333 days for 
Jupiter’s period, and amplitude 0%25. No cur- 
rently known source of systematic error would 
account for such a period. The residuals for the 
other outer planets contain no period as marked 
as the 4520 day period in Jupiter. A secular 
increase in the eccentricity of 1/3 X 10-® per 
century, together with a motion of the perihelion 


Saturn does give some indication of a 33 year 
period. A 400 year period in Neptune would pre- 
sumably be indistinguishable from a constant 
displacement and so would have been removed 


e— = — cos d@—o)+ 2 sin d@—w). (9) 


. constant force producing an acceleration with 
omponents F, and /; along the major and minor 
xes respectively then leads to secular variations 


de ar, 3 fr amounting to 5’’/century, will, if coupled with a 

. dt 20 term whose frequency is equal to Jupiter’s mean 
nd (10) motion, produce such a period. If such secular 
5 da _ eae perturbations were produced by a constant anom- 

dt DD alous acceleration of Jupiter, the magnitude of 


the acceleration required is'about 10~® times the 
centripetal acceleration of Jupiter. The secular 
change in eccentricity is too small to be statis- 
tically significant. The data for:the other outer 
planets are inadequate either to confirm or deny 
the existence of such an acceleration in their 
motion. 
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For Jupiter a force along the major axis and 
roducing an acceleration whose magnitude is 
bout 10~® times the centripetal acceleration vn 
ould produce the required motion of the peri- 
elion. Such a force acting on the other outer 
lanets would produce larger amplitudes propor- 
onal to 1/v and longer periods. The periods and 
mplitudes to be expected are shown in Table II. 


TABLE Ii. PERTURBATIONS PRODUCED BY A CONSTANT 
ANOMALOUS ACCELERATION 


Planet Amplitude Period of this work. 
yupiter 0%25 12.4 years REFERENCES 
Saturn 0.34 33 
Uranus 0.48 120 Clemence, G. M. 1951, M. N. 111, 219. 
Neptune 0.60 400 . 1954, Astr. Papers Amer. Ephemeris 13,367. 


Clemence, G. M. and Brouwer, Dirk. 1955, A. J. 60, 118. 


The residuals for Saturn, Uranus and Neptune 
either confirm nor deny the existence of such 
erturbations, although the auto-covariance of 


Eckert, W. J., Brouwer, Dirk and Clemence, G. M. 1951, 
Astr. Papers Amer. Ephemeris 12. 

Hertz, H. G. 1953, Astr. Papers Amer. Ephemeris 15, 169. 

Tisserand, F. 1889, Traité de Mécanique Céleste, Vol. I. 


STABILITY OF STRAIGHT-LINE SOLUTIONS IN THE RESTRICTED 
PROBLEM OF THREE BODIES 
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Abstract. About thirty orbits of infinitesimal particles near Lz and L; are computed for » = 0.1 with the IBM 701 com- 


iter of the University of California. 


The results do not favor accumulation of material near these points of equilibrium. 


table of newly computed coordinates of Li, L2 and L; for various values of u is given. 


1. Introduction. In the study of close binary 
‘ars the problem of concentration of material at 


Lagrangian points L, and L; is of considerable 
interest. According to the first order theory, one 
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would expect stable orbits around L»2 and Ls; 
under some special initial conditions. But nu- 
merical calculations are needed for estimating 
the effect of the second and higher order terms 
in disturbing the stability of these solutions, and 
also to delimit the region of tolerance for the 
initial conditions under which the orbits are 
stable. If the orbits are sufficiently stable and if 
the initial conditions are not too severe, only 
then will a significant amount of material accu- 
mulate near the equilibrium points in question. 
The IBM 7o1 digital computer is ideal for these 
computations. Gould (1957) had written a pro- 
gram for the machine for calculating orbits of 
particles ejected by prominence action from the 
surface of one component of a close binary sys- 


tem. Her program is used here for obtaining som 
representative results for the case of w = 0.1. 

After these computations were made, h 
author’s attention was drawn to a rigorous proo! 
based on Lyapunov’s method, given by Duboshi 
(1938) for the general nonstability of the periodi 
orbits near the points Ly and L;. It is hoped tha 
the numerical results presented in this paper wil 
serve as an illustration of the actual conse 
quences of nonstability for varying startin 
conditions. : 

2. Initial conditions. The stability of straigh 
line solutions in the restricted problem of tht 
bodies is discussed by Moulton (1947). His’ 
sults are collected in this section for clarifyi 
the procedure of the present computations. F« 
lowing Moulton’s notations, the straight-lin 


TABLE I. CONSTANTS OF MOTION AROUND Lz AND Lz 


va 


Point Aj o/c —co e fa 
Le 2.5067 42371 0.63032 14206 —4.3904 34977 3-777 0.9254 
L3 1.0916 91979 0.53412 85365 —2.1716 66505 5.834 0.8684 
* The period of revolution of the system is 27 in the same units. 
TABLE II. FIRST ORDER ‘EQUILIBRIUM’ ORBITS AROUND L2 AND L3 
Initial conditions 
Point Orbit Xo! =X — XL Mo" =o iq! =Xq do’ =o 
L2 Liner -++0.1000 +0.1000 +0.0630 32136 —0.4390 43024 « 
Il 2 —0.0500 —9.0500 —0.0315 16068 +0.2195 21512 7 
Il3 +0.0100 +0.0100 +0.0063 03214 —0.0439 04302 i 
Il 4 +0.0010 —0.0020 —0.0012 60643 —0.0043, 90430 i 
Il5 +0.0001 +0,0002 +0.000I 26064 —0.0004. 39043 
L; Ili 1 —0.1000 0.1000 —0.0534 12854 +0.2171 66650 3 
III 2 —0.0100 —0.0200 —0.0106 82571 +0.0217 16665 
Ill 3 —0.0010- —0.0020 —0.0010 68257 +0.0021 71667 
Ill 4 —0.0001 — 0.0002 —0.0001 06826 +0.0002 17167 
TABLE III. ‘NONEQUILIBRIUM’ ORBITS AROUND Lz AND L3 
Initial conditions 
Point Orbit Xo! v0! Ai’ Ado! 
Le II 31a +0.0100 +0.0100 +0.010000 +0.010000 
II 31b +0.0100 +-0.0100 —0.010000 —0.010000 
II 32a +0.0100 +0.0100 +0.001000 +0.001000 
II 32b +0.0100 +0.0100 —0.001000 —0.001000 
II 33a +0.0100 +0.0100 +0.000100 +0.000100 
II 33b +-0.0100 +0.0100 —0.000100 —0.000100 
II 51a +0.0001 +0.0002 -+0.000100 +0.000100 
II 51b +0.0001 +0.0002 —0.000100 —0.000100 
II 52a +0.0001 +0.0002 +0.000010 -+0.000010 
II 52b +0.0001 +0.0002 —0.000010 —0.000010 
II 53a +0.0001 -+0.0002 +0.000001 +0.000001 
II 53b +0.0001 +0.0002 —0.000001 — 0.000001 
L; III 31a —0.0010 —0.0020 +0.010000 +0.010000 
III 31b —0.0010 —0.0020 —0.010000 —0.010000 
III 32a —0.0010 —0.0020 +0.001000 +0.001000 
III 32b —0.0010 —0.0020 —0.001000 —0.001000 
III 33a —0.0010 —0.0020 -++0.000100 +0.000100 
III 33b —0.0010 —0.0020 — 0.000100 —0.000100 
III 34a —0.0010 —0.0020 +0.000010 +0.000010 
III 34b —0.0010 —0.0020 —0.000010 —0.000010 
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Figure 1. Orbits of infinitesimal particles around Lz and Ls; for u = 0.1. The xy system of co- 
ordinates is rotating anticlockwise with unit angular velocity in a fixed frame of reference. Num- 
bers with marks on the curves denote units of time elapsed since the particle was started. 


lutions occur at x = Xo:, y = 0, 2 = O, where 
; are the roots of equation: 


28 a Bef 
7 Steers a 
x — Xo 


heat 


ae coordinates are referred to the rotating frame 
reference and the units are canonical. If 
= x9; + x’, y = y’, 2 = 2 are the coordinates 
an infinitesimal particle placed close to any 
ie of the three Lagrangian points Li, Le and Ls, 
that x’, y’, 2’ are small, the equations of motion 
the particle are reduced to: 


ax! dy’ ; 

ee de (1 + 2A,)x 

dy’ dx’ 

Se bear = (1 — 2A;)y"?, (2) 
dz! : 

nary” 


where 


: Saal! 2 


~ [Gos — x1)?]}? zy E (xo; — x2)? |? (3) 


A; 


The equation in the 2’ coordinate is independ- 
ent of the others; it gives periodic motion in the 
z’ coordinate with a period of 27/A,;. Limiting 
ourselves to the x’y’ plane, a nontrivial solution 
of the form x’ = K e, y’ = L e™, is possible if 


M+ (@-—4)V+U04 4; - 2A?) =0. 


Since (1 + A; — 24,7) <0 for all » < 0.5, we 
have 
M=-+P, or =— Q, 
where P and Q are positive real numbers. There- 
fore 
M=—Aw=+ VP, and A=— DM _ 
=+ iV0 =+ ic. 


} 
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+130 + 


Figure 2. Orbits of infinitesimal particles around L2 for » = 0.1. The whole units of time are marked 


About thirty orbits around Ly and L; are 
puted. Nine of them were started at differ 
distances from Le, or Ls, with velocities de 
mined by Eq. (5), which would make the or 
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+0.08 

+0.04 

0.00 

-0.04 

-008 

F115 +1.20 41.25 
on the curves, except where fractions are indicated. 
Also, 
Ho (ea ule 
Ly= ue i Kj = 6; K;; 
J 

and again c¢; = — Co, ¢3 = — c4'= 7 cy where 


C= (co? +1 + 2A,;)/20. 


In order to get periodic solutions, we set Ky 
Ke = 0 and, therefore, also L; = Ly = 0. Then 
the initial conditions become: 


dxo'/dt = (a/c) yo', dy /dt (5) 


Under these conditions the particle presumably 
follows a periodic elliptic orbit centered on the 
Lagrangian point in question. It would have a 
period = 27/o and e = V(2 — 1)/c?. In these 
calculations we have taken » = 0.1, R, = 0.1, 
Ry, = 0.0125, Ri and R: being the radii of the 
two components. Table I gives A, o/c, —oc, P 
and e for Lyand L;3. 


(4) 


a Coxo 


stable according to the first-order theory. Th 
initial positions and velocities are given 
Table II; the units are canonical. The rest of 
orbits were started at x’ = y’ =+ 0.01 % 
x’ =+ 0.0001, y’ =+ 0.0002, for Le; 


’=— 0.001, y’ 0.002, for Ls, with com 


xX = 
ponents of velocity differing from the velocit 
of the first-order theory by amounts Az’, M& 
given in Table III. Figures 1 to 4 represent # 
results: Figure 1 shows the orbits on a sm 
scale, Figures 2 and 3 represent the orbits arou 


L, on larger scales, and Figure 4 is a large s 


- drawing of the orbits around L;. Numbers v 


the marks on the curves denote units of tin 
elapsed since the particle was started. 


car: 
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3. Results. Our conclusions are: 

(i) Over a wide range of initial distance from 
», even with the proper initial conditions, the 
urticles make only one revolution around Lz, 
}fore they go off. Particles moving near Ls; 
jake two rounds if the distance from Lz, is 
\fficiently small. The orbits are partly elliptical 
ad their periods and eccentricities are in fair 
sreement with the predictions of the first-order 
ear. 

(ii) The incipient stability is destroyed even 
‘the initial velocity components differ by about 
te per cent of the velocity components given by 
a. (5). For bigger deviations the orbits depart 
ore and more quickly from the elliptical orbits 
the first-order theory. 

(iii) These results are not favorable for ac- 
(mulation of material near the Lagrangian 
vints L, and L3. Gould (1959) had found that 


02 


22 
Sy, 
: v3 
Dol 73 3 
e 
| Ss 
| vs 
5 
woo!_. X-axis 
| 
[oXe) 
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the particles leaving either component of a binary 
system stay near L» and L; only for a short time 
during which their orbits are changed from di- 
rect to retrograde motion. 

(iv) Most of the computed orbits are fun- 
neled into orbits close to each other, which sur- 
round both components of the binary system. 
This result gives support to Gould’s conclusions 
regarding the formation of rings in close binary 
systems. 

(v) A few orbits were computed for the case: 
= 0.5. The results were similar to those for Lz 
in the case of uw = 0.1. 

4. Coordinates of Li, Le and L3. In computing 
A;, o, and c, the coordinates of L, were taken 
from the work of Kuiper and Johnson (1956). In 
order to verify their results and also for comput- 
ing the position of L3 for » = 0.1, a small pro- 
gram was written for obtaining x; for the three 


35 


49 245 


+1.262 


+1.261 


+1260 


Figure 3. Orbits of infinitesimal particles around L»2 for » = 0.1. The whole units of time are marked 


on the curves, except where fractions are indicated. 
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-1.050 -1045 -1040 O35 max 
Figure 4. Orbits of infinitesimal particles around L; for 4 = 0.1. The whole units of time are 
marked on the curves, except where fractions are indicated. 
straight-line solutions. This program was used :5 — (3 — nu) pit + (3 — 2u) pb 
for computing the positions of the three straight- — (er — 291 + 1) =m 
line solutions for other values of yu. The results { 
are given in Table IV; they are good to within a From this we arrive at the iteration formula 
few units in the last place. j 
If p: is the distance of L; from yp, we obtain — u(I — ps)? | G 
from Eq. (1), PY" LB = 2H) — pif — #) — pil | 


oe 


).0000 OOIO 
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TABLE IV. COORDINATES OF Lj, Lz AND L3 FOR VARIOUS ws 


| KB m2/m1 XL, 
).5000 0000 1.0 -+-0.0000 00000 
».4444 4444 4/5 0.0785 05894 
).4000 0000 2/3 0.1416 17526 
1.3750 0000 3/5 0.1773 42630 
).3000 0000 3/7 +0.2861 29783 
).2857 1429 2/5 +0.3072 33226 
').2500 0000 1/3 0.3607 43428 
).2307 6923 3/10 0.3900 97486 
»).2000 0000 1/4 0.4380 75959 
1500 0000 3/17 +0.5197 40359 
).1000 0000 1/9 +0.6090 35110 
).0500 0000 1/19 0.7152 25350 
).0200 0000 1/49 0.8034 65629 
),0100 0000 1/99 0.8480 78713 
).0030 0000 +0.9003 73655 
).0010 0000 +0.9312 86976 
-),0003. 0000 0.9540 10057 
).000I 0000 0.9680 65213 
).0000 3000 0.9785 81345 
1.0000 1000 +0.9851 26720 
).0000 0300 +0.9900 30459 
).0000 O100 0.9930 81618 
).0000 0030 0.9953 65315 
+0.9967 85658 
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XL XL 

+1.1984 06144 —1.1984 06144 
1.2170 25950 1.1785 83395 
1.2308 13769 1.1620 45267 
1.2380 38366 1.1525 24009 
+1.2567 34696 —1.1232 05596 
+1.2596 66516 —1.1175 14574 
1.2658 58102 1.1031 66848 
1.2684 02205 1.0953 78413 
1.2710 48690 1.0828 39465 
+1.2703 34073 —1.0622 98631 
+1.2596 99833 —1.0416 08909 
1.2280 93667 1.0208 26334 
1.1800 77904 1.0083 32894 
1.1467 65042 1.0041 66612 
+1.1002 85368 —I.001I2 49999 
+1.0699 16097 —1.0004 16667 
1.0468 25818 I.000I 25000 
1.0324 25189 1.0000 41667 
1.0216 68072 1.0000 12501 
+1.0150 02015 —1.0000 04168 
+1.0100 30190 —1.0000 OI251 
1.0069 48468 1.0000 00417 
1.0046 48045 1.0000 00125 
+1.0032 21637 —1.0000 00042 


Then x1, =1—yp— p:. Similarly for Ly we 
ybtain ps, the distance from uw by the iteration 
‘ormula: 


| u(L + pe)? ( 
| tee earcemess a 
ied XL, = I — p+ pe. 


If ps is the distance of L3 from (1 — u), we 
| p3 by substituting (1 — yw) for » in Eq. 
(7). Thus, 

| iigaecf) (Li ps)? | (8) 
i? | © 2n) + af (2 +) + ps} J’ 


1 


with Ee eh —. Pe 
_ The results given in Table IV are obtained 
from Eqs. (6), (7) and (8) by iteration from a 


starting value of p; = 0. They are in agreement 
with the figures given by Kuiper and Johnson 
(1956), and Kopal (1954), up to the number of 
places given by them. There is a discrepancy of 
58 in the last place of xz, and x1, for » = 2/7, 
(m2/m, = 2/5), given by Kopal. 

I am grateful to Prcfessor O. Struve for sug- 
gesting this problem and to Mrs. Nancy L. 
Gould for the use of her program. 
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RELATIVE PROPER MOTIONS OF STARS IN THE REGION 
OF THE OPEN CLUSTER NGC 2281 * 


By S. VASILEVSKIS anp ALBERT G. A. BALZ, JR. 
Lick Observatory, University of California, Mount Hamilton, Calif. 
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Abstract. Relative proper motions of 127 stars in the region of NGC 2281 were determined from two pairs of plates 
taken at the Allegheny Observatory. A difference in epochs of over 40 years gave annual proper motions with p.e. +0*00¢ 
in each coordinate, after the measurements had been corrected for a slight magnitude effect. 4 i 

The total membership in the field of approximately 50’ X 40’ down to a plate limit of about photographic magnitude 
13 is approximately 56; the probability of membership was derived for individual stars. 


The open cluster NGC 2281 has been classified 
by R. J. Trumpler (1930) as Type I3p in appear- 
ance and Ia in spectral composition. As far as the 
authors are aware, no astrometric investigation 
has been done in the region of this cluster. The 
present study is an attempt to separate the mem- 
bers of the cluster NGC 2281 from the field stars. 

Plates and measurement. Two  first-epoch 
photographs taken by Trumpler with the 30-inch 
Thaw refractor at the Allegheny Observatory 
were repeated with the same telescope, and both 
plate-pairs, listed in Table I, were kindly made 
available to us by Dr. N. E. Wagman, Director. 

Each 8 X Io-inch plate was measured in direct 

and reverse positions, with two bisections on each 
measurable star image in each position. A total 
of 127 stars was measured; they are shown in 
Figure I. 
’ Plate-pair 1 has three exposures of different 
length; all images of the longest exposure were 
measured, and brighter stars were measured also 
on shorter exposures, with a sufficient number of 
stars common to two consecutive exposures. The 
measurements of these common stars were used 
to compute constants of linear reduction from 
systems of the shorter exposures to that of the 
longest one. ae 

Although plate-pair 2 contains more images 
than any exposure on plate-pair I, only stars 
common to both plate-pairs were measured, and 
two is the minimum number of measurements 
for each star. The number of measurements for 
individual stars is given in column 9 of Table II. 


TABLE I. PLATES MEASURED 


. Plate 
Pair No. 


Date 
I 4350 1915 Dec. 22 
89216 1957 Mar. 10 
2 4262 1915 Nov. 30 
89063 1957 Feb. 2 


* Lick Observatory Bulletin, No. 561. 
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All the plates were measured and most of re 
ductions were carried out by one of us (Balz)._ 

Reduction of measurements. Differences 
measured rectangular coordinates on plates of 
both epochs were formed separately for eac 
exposure-pair. On plate-pair 1 the differences fo 
shorter exposures were reduced to the system of 
the longest exposure, using measurements of 
stars common to two consecutive exposures 
Then the mean for all the exposures was taken 

The frequency distribution of the difference 
in x and y was examined for the purpose of 
lecting the reference stars from possible - 
stars with small proper motions. Twenty-five 
stars with nearly uniform distribution over h 
plate were selected, and they are marked wi 
asterisks in column 1, of Table II. As found lat 
most of the reference stars were probable met 
bers of the cluster; a possible systematic effect 
of this sélection was, therefore, taken into até 
count, as explained. below. a 

Vector-points of proper motions showed a p 
nounced concentration, with an indication of ¢ 
magnitude effect. The vector-points of the 
brightest and the faintest stars seemed to be 
placed in the same direction, with reference t 
the points of the medium bright stars. This typ 
of error can be explained as a guiding error. 

Let us assume that for a certain fraction 71 ¢ 
the total exposure time the image was off-set in| 
certain direction from the guided position, ant 
for another fraction 7, it was displaced in th 
opposite direction. 


=A 


If +; > 72, and each of th 


Exposures Interval 


(sec.) Seeing (years) 
5, 25, 150 4,3 41.21 
ype Sora Co)s) AA. 

600 Rh RS 41.18 
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Figure 1. Identification chart of stars measured. 


actions is small compared with the total ex- 
psure, the exposures of the displaced images will 
Il below the emulsion threshold for the faintest 
ars. With increasing brightness 71 will become 
‘ng enough to have an effect in position, but the 
{posure 72 will still have no influence until the 
crease of the brightness is sufficiently large. 
hen the exposure 7; will start to act in the op- 
site direction compared with 71, i.e., it will 
ove the centroid of the image towards the 
‘stem of the faintest stars. This type of magni- 
ide effect was determined and taken into ac- 
unt as follows: 

|A circle with radius of 0730, approximately 
vice the mean error of a centennial proper mo- 
on, was drawn around the condensation of the 
sctor-points, and all the stars with the points 
ithin this circle were assumed to be probable 
embers of the cluster. For each of these stars 
1 equation was written, separately for x and y, 
the form: 


100u, = 100u’ + Aa + Bb + Dd+ Ee (1) 


where 100,, is the centennial relative proper mo- 
tion of the cluster and 100u’ that of a suspected 
cluster member. Other symbols denote: A = x/ 
100, B= y/100, D=10 (s — 0.200) and 
E = D?, where x and y are the approximate co- 
ordinates of a star, and s is the diameter of the 
image of the star on one of the photographs; 
x, y and s are in millimeters. The constants d and 
é are introduced to account for the suspected 
magnitude equation, as described above. Sepa- 
rate least-squares solutions of (1) in x and y gave 
the following values with corresponding prob- 
able errors: 


+o01I = 100p’,—0.05A +0.01B+0.01D+0.08E 
+0.02 0.04 0.05 =0.02 =£0.04 


—0"15 =100p’,—0.01A —0.12B+0.04D —0.10E 
=-0502 =-0:04 =50.05 =£0:02 =-0:04 
(1’) 
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TABLE I. RELATIVE PROPER MOTIONS 
No. C Me x y 100, 1004, Py n No. C mi, x y 100, 1004, Py 
(1) (2) (3) (4) (5) (6) (7) (8) ~— (9) (1) (2) (3) (4) (5) (6) (7) (8) 
1 12.7 -25.4 415.4 -006 -0:19 0.95 2 41 13.2 - 3.8 +12.6 41.09 -0.05 
2 11.5 23.6 +11.7 + .84 - .80 3 42 61° 13.0 3.3 =12:0)\+ o12er noe 
3 10:1 23.2 414.3 4+ .62 -2.06 3 43. 73° 11-0. 3.10 = 3.5. = 910 meenoteeeom 
4 13.2 23:2 +10)2 4.32 4 558 2 44* 74 10.7 3.0 = 2.0) = s0Gu-Odmuenon 
5* 12-5 92)5) = OSe= KObt eNOL Olen 45. 76 ° 9.7 ~/3.0 =00,1 4 201 cecOmeege 
6 12.4 -21.6 +1.8 - .50 + .20 2 46112. 11.8 = 2.7 =19.'7 1G ae 
7* 1125" 10.7 = 8.5 4.46 =), 06m403) 2 47 12.9 2.6 414.6" 4130 36 
8* 10:8 18.8 415.9) 4.5119) 324) 9294) US 48° 78 411.8 2:1 + 2.8 4.16 seus OGumenan 
9 10.9 17.9 -10.5 -1.59 +2.34 3 49 23 10:3 1:4 - 6:4 + $14 ptomeege 
10 11.9 LS er mOn Se beerr ete 2 50. 131 - 11:7.~.1.2) + 354-4) Soe 
11 Tout 1720) -=1 2940-8 Oe o5 2 51 82 fist. - 1.1 + 1.9) 412m oeae 
12 19.5 1619) =10,4) + 201 42501 2 52-59 10.7 - 0.8 -11.5 +.05 -.14 .96 
13 1100) 21688) ONG ets 22 1G 3 68-26 10.8 . 0:0! 5= 3./0) — aetna 
14* 1950). 16) 0p 14:9 Ga coe oie 54 80 --10.5 0.0. 0:0) + O7meuOMneg 
15 12.8 16.7 -14.3 + .36 +.75 2 55 20 ° 9.3 + 0.1 - 6.3 4.05 mounaamGG 
16 13020215299 =n2- Oni=t- 04a n sta Bs ee 56. 24 12.2 + 0/2 = 479) Sivoomenere 
17 11.7 14.9 + 5.6 4.25 -5.69 3 57 5° 28. 10.2. 0,4<-92. 9) 2duamnG mer 
18 9.0 4407 4 Oot) 4-02 en 08 BseTOe 58 21 9.2. "0/5. ='5.70 4) 417 jae Oa 
19 A1esveisuds 2s 007'=.09 96S 59 «19 «98. "0.8 = 6,3. = 204 
20* 10a. 1249.4) 720524 2 Sou One olaLS 60. 14. -18.0-= 0.9. - 852) +659 ies Omens 
21 1258) =1259 s=1004 + 42m 02s 1One2 61° 18 -10.4-+:1.0 = 6)4 4 11 eel Oeeoe 
22* i128) -112.9° 41859 = 104 c=5 2on eos) ve 62-17 «19.0: -.1.2) =.5)6) 412 ee o eon 
23 hey sb eclecs Kioeey Su) ech a 63. 134. 8.6 1.3. + 7.4 = 410 = ouenon 
24 12393). 129) C1907 a Ghee be 2 64-17 10:8 =~ 1.3 = 5.5 +..08 =i Saumeag 
25 1126.) 1201) 812.0) Ee Soe eAD 2 65 55 $8.8 | 1.5 +1269 (4.02 aomeOn 
26 13.4 -10.9 -18.1 +.02 -1.46 2 66 58 11.8 +1.6 -10.4 + .24 -.16 .83 
27* 1057 10:8) 20:5) =4.055— 401") 9311.9 67* 141 12.2 1.6 -20.0 +.63 - .08 
28 Ty CES Sy so Sivas Gh) Sl 8 68* 9.8 2.1 413.0" ='-07) = Ocmeeae 
29% {2s4= 8 5tats oe —e tOk eae 2 69 57 12:7 2.1 -1053 4.32) Oomeoe 
30° 1192 (9:1 8:0n-1619) =7.09 = 06, 94 4 70 (56 12:0: -2:1° =1050, =c65memeet 
31* 1055) c= o7 Te 16,5042 200 0S odd 7. 12 11.3 + 2.3 - 9.2 +.75 -3.60 
$2. 417 1244 6) 110; 6 ps3 90 2, 72 11.1 2.4 419.0 41.10 + .56 
33 12.9 6.4 414.8 +2.81 41.25 2 13* 86 9 12.7) 23554 4:2) 2). 1S eee 
S470 1258) Gein — 6) be poo ak 2 14°41 10:0. ° 2.6 —15:8 0 y2medes 
SO al 25nd G6 obi oa aoe TT 3 75 140. 11.3 2:6 -18)4 = 912semos 
B60 104° 912500 "5.4; 1940, = 18k 20 80m 2o 76 88 12.8 +3.1 +5.7 +.50 -1.69 | 
37 12030 405 io). idee 04 pansy oun 7] 52 11.4 3.2 -13.1 + .26 -.06 75/0 
38 115 9.8 4.4 -12:% = .03 =.18 96 °4 78 4 9.9 3.3 - 4/2 -.26 231) Sei 
395) 6550 1256, 4996.65 = 03a A Osmesnice 79 5S 9.7 3.9 = 5.2 =.05 —- 143 0;¢QumE 
40) 1755) 2.958) = S.9umwteon en oe 710940 a4 80 89 11.4 4.1 +5.3 +.22 +.35 


The solution (1’) shows by the size of the co- 
efficient of E that the magnitude effect described 
above is probably real: Systematic effects in 
scale and orientation, caused by the magnitude 
effect, or by selection of the reference stars, are 
also taken into account by (1’). The corrected 
relative proper motions are plotted in Figure 2 
and they are listed in Table II in columns (6) 
and (7); the concentration of vector-points after 
the correction is better than before. The mean 
error of a centennial proper motion of an average 
weight is o =+0"13 in each coordinate; the 
corresponding probable error is +0"09. 


‘done in the earlier paper. It has been assume 


i 


Cluster membership. The probability for clus 
membership has been derived for individual st 
in a way simplified as compared with an earl 
paper (Vasilevskis, Klemola and Preston 1958). 
Figure 2 shows that, although proper motions ¢ 
the field stars are not distributed uniform 
around the condensation of points, the num 
of field stars seems to be too small and the scatte 
in proper motions too large, for an elaborate der 
vation of their frequency distribution, as 


here that the vector-points of the field stars ar 
distributed uniformly within a circle of radius 4 
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TABLE II. (continued) 
x y 100p., ee p n No, C¢ mE x y 100, ECO. p n 
(4) (5) (6) (7) (8) ~— (9) (1) (2) (3) (4) (5) (6) (7) (8) ~— (9) 
+ 4)1 - 7.3 40/03 -0.11 0.96 2 106 13.0 +15.0 -11.5 -0!'26 +0.39 2 
4.3 - 4.3 - .05 - .25 .94 4 107 L229 pep Onl Sele ON 4+ 0.43 2 
4.8 - 3.3 +.13 - .34 .84 4 108* LOSO L516 4 7227 +2536 4 07 0:14 — 3 
4.9 -18.9 - .31 - .50 .03 2 109* 12.1, 15.7 - 9.8 00 - .07 9523 
4.9 - 3.2 +.21 - .09 .87 3 110* 12.6 15.9 - 4.4 + .09 - .22 95) °2) 
+5.0 - 6.7 +.01 - .16 96 4 111 12.6 +16.7 +11.2 - .62 +1.34 2 
Dente = 12 .< 207 49 2 112 9.7 16.8 +19.8 - .67 + .26 4 
‘| 5.4 -20.2 + .50 + .46 3 113* ee Lili Ou hme 02 + .08 83 2 
| 5.4 -11.4 - .07 + .14 63. 2 114 12.8 18.2 +11.8 + .08 +1.01 2 
| 5.4 +22.2 + .68 - .03 2 115* 12.5 19.4 -16.9 - .43 + .01 04 2 
1 6 11.2 +5.5 - 5.5 Nooo =970k 232 2 116 10.3 +19.4 - 8.9 -1.35 - .85 4 
2 106 10,2 6.0 -13.9 - .08 - .15 .95 3 117 878. 192-8 =10257 —2507 =-..16 95 4 
3 41 le} 7.7 - 6.8 - .04 + .11 {ieee} 118 13.0 20:3 - 3.8 + .85 + .81 2 
4 38 12.9 8.9 - 3.5 +.05 + .88 2 119 12.4 20.8 - 6.5 + .95 -3.02 2 
mt 104 10.1 9.5 -10.5 + .16 - .02 .88 4 120 11.0 22.5 +12.1 + .07 -3.80 3 
12.1 +9.6 + 9.8 - .42 -2.03 2 121 12.9 +23.5 + 8.5 + .43 - .28 .08 2 
ng 10.5 11.0 -18.2 + .80 - .59 3 122 12.4° 24.8 -11.2 + .03 + .83 2 
(8 137 12.9 11.2 -15.2 + .03 + .76 2 123* i255 2620) 497.8" = 225 + .24- 9320) 33 
‘9 10.1 11.5 +14.0 + .87 -1.35 4 124* 11.0 26.6 +16.4 + .21 - .44 .09 3 
0 11.2 11.7 +16.7 - .23 + .44 3 125* 1152 26.7 = 2.9 - .10 -.08 .94 3 
1 94 13.2 411.8 + 2.4 - .86 + .11 2 126 12.0 +29.0 - 1.3 + .63 - .51 3 
2* 12.3 12.6 +19.8 - .45 - .13 .05 2 127 1258 29.1" = 953 "+ 215 +-.03' (0.84 22 
3 103 11.5 13.0 - 9.9 - .99 -2.20 3 3 
4* 12.5 13.6 415.7 + .44 - .36 0.02 2 
5 97 13.0 14.3 - 5.2 + .32 + 82 2 
| EXPLANATION OF COLUMN HEADINGS 
| Col. (1) : star numbers; asterisks denote reference stars. 
| Col. (2) : star numbers by Cuffey (1938). 
Col. (3) : photographic magnitudes. 


with star No. 54 as origin. 
\ (6) and (7) : centennial proper motions. 


(9) : number of both-epoch measurements. 


found the center of the condensation, where 
‘/=+ 0"13 is the mean error of a centennial 
soper motion. The center is defined by the solu- 
on (1’), and is marked by a cross in Figure 2. 
|Under the assumption that the vector-points 
' the cluster members represent a_ bivariate 
»ormal distribution, the fraction of the total 
)pulation within a circle of radius 7, expressed 
| units of the standard deviation, o, can be 
jund from 


= 
y= uF Ree de = 1 — ere (2) 
0) 


addition to the circle of r = 4a in Figure 2, 
iree smaller circles of r = c, 20 and 30 were 
rawn but not shown in Figure 2. Then Ay, the 
action of the total population within the small- 


Col. (4) and (5) : approximate rectangular coordinates in minutes of arc, oriented for 1950.0 and 


(8) ; probability of a star being a cluster member. 


est circle and within the zones between each two 
consecutive circles, was computed and tabulated 
in column 3 of Table III. The circle r = 4o in- 
cludes all the members of the cluster. 

Column 4 of Table III gives the fraction, Ak, 
of uniformly distributed vector-points of field 
stars within corresponding zones, and column 5 
shows the number of actually counted vector- 
points in each zone. If JN is the total membership 
of the cluster and 7 is the number of field stars 
with vector-points within 7 = 4c, then each zone 
gives an equation 


AvN + Akn = S, (3) 


and the corresponding least-squares solution 
yields N = 56.3 and m = 18.3. These values 
inserted into (3) give columns 6, 7 and 8 of 
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TABLE II. DISTRIBUTION OF VECTOR—POINTS 
r y Ay Ak S AvN * Akn & 
(1) (2) (3) (4) (5) (6) (7) (8) 
o +394 394 .062 25 22.2 Doe 23 
20 .865 471 .188 29 26.5 3.4 30 
30 .989 .124 a2 ma 7.0 Eye 13 
4o 1.000 OIL 438 10 6 8.0 9 
1.000 1.000 75 56.3 18.2 75 
Table III. A comparison of columns 8 and 5 cluster member, was computed for narrow zon 


shows a good agreement and justifies our as- 
sumptions on the distribution of vector-points. 
The ratio of the expected number of vector- 
points for cluster members to that for all stars, 
i.e. the probability of a point belonging to a 


Figure 2. Relative proper motions of stars measured. 
Proper motions of 13 field stars were too largé to be plotted 
within the‘limits of the figure. 


and plotted graphically. Then the probabil 
b:, were read from this graph for each vecto 
point within Sh AG and listed in Table I 
column 8. ~ 

Photographic ue Mpg, Were delond 
only approximately. Magnitudes by Cu 
(1938) were used as standards and _ plotte 
against the diameters of images on one of tf 
plates measured. Then the magnitudes for a 
the stars were read from this calibration gr 
The probability frequency distribution of clust 
membership in m,, is shown in Figure 3. 


M 
=D 
° 


- wn ww am ODN @ 


° 


9 10 iT 12 13 


Figure 3. Frequency distribution of a 
cluster members in mp9. if 
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SOUTHERN HEMISPHERE PHOTOMETRY VI 


THE COLOR-MAGNITUDE DIAGRAM OF NGC 458 AND THE ADJOINING REGION 
OF THE SMALL MAGELLANIC CLOUD * 
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| 
| Abstract. The first color-magnitude diagram of one of the so-called ‘‘blue globular clusters’ in the Magellanic Clouds 
|sexamined. Although the luminosities of the stars it contains are very much like those of the Pleiades in our own Galaxy 
jie color-magnitude diagram is very different. There are a large number of yellow giants, a narrow and well-defined 
lertzsprung gap, and a discontinuity of about one magnitude between the main sequence and those giants. 

Because of these striking differences between this Small Magellanic Cloud cluster and clusters of approximately the 
me age in our own Galaxy, it is concluded that the chemical composition of the stars in NGC 458 is different. 

| The color-magnitude diagram of the field adjoining NGC 458 is quite similar to those of the fields shown in earlier 
apers of this series for two other parts of the SMC. Now it is possible, however, to identify the major kinds of stars in 
Il these fields to be primarily a mixture of the older stars found in the clusters like NGC 419 and NGC 361 and those 


und in the younger clusters like NGC 458. 


1. Introduction. In 1951, Gascoigne and Kron 
1952) found that the integrated color indices of 
‘I star clusters in the Magellanic Clouds fell 
bout equally into a red group and a blue group. 
“he red group was tentatively assumed to corre- 
pond to globular clusters in our own Galaxy, 
ind representatives of that group have been ex- 
‘mined in earlier papers of this series: SHP III 
NGC 419; Arp 1958c) and SHP V (NGC 361; 
\rp 1958d). c 

The blue group was even more interesting be- 
ause, though it was assumed that these clusters 
yere composed of blue stars such as occur in 
‘pen clusters in our own Galaxy, their geomet- 
ical structure more or less resembled that of the 
lobular clusters (see Gascoigne 1954, where 
VGC 458 and NGC 1866 are listed as globular 
llusters). Their appearance is typified by NGC 
58, shown in Plates I through III in the present 
yaper, or by NGC 1866, an even richer object 
h the LMC. 
| Since NGC 458 was the brightest of these 
ibjects in the SMC and favorably situated near 
he northeastern edge, it was chosen as one of 
he three areas in the SMC where faint calibra- 
ion stars would be measured. The analysis of 
his region completes the report on the three 
3MC probes to faint limiting magnitudes. Only 
inalyses to brighter limiting magnitudes in the 
nore central NGC 330 and NGC 371 regions will 
»e reported in forthcoming papers of this series. 
2. Magnitude standards. Exactly the same tech- 
liques were used in the present NGC 458 area 
0 measure the photoelectric calibration stars as 
vave been described in SHP III and V. The 
yhotoelectric values listed in Table I and identi- 
ied in Plate I represent an independent sequence 
‘onnected to the usual network of local standards 


reported in SHP II (Arp 1958b), but the ob- 
serving priority assigned less observational weight 
in NGC 458 than to the primary NGC 419 
sequence. The adopted values in the last two 
columns of Table I were obtained from the mean 


TABLE I. MAGNITUDE SEQUENCE IN NGC 458 


Photoelectric measures Adopted values 


Star V B 4 B 

a 9-74 10.16 9.74 10. 16* 
d 9.79 10.80 9.79 10.80* 
b 10.87 Tin33 10.87 Tin35) 
G 11.53 12.06 11.53 12.07" 
e 12.54 12.91 12.47 12.91 
g 12.58 13.25 12.58 13.24 
a5 14.37 14.38 14.38 14.10 
br 15.21 15.06 15.19 15.12 
el 16.66 16.37 16.33 16.37 
b2 16.34 17.98 16.50 18.10 
c5 16.72 17.45 16.66 17.56 
ZN Ss oR See tone 17.59 16.73 17.49 
aq OOS tees 16.99 18.347 
GLO abate see 18.31 T7231 L7eOM se 
c3 17.46 18.68 17.55 18.56 
DME Nd pe Sth edie dren 0 17.80 18.98 
CL URS) Sl Of em tae ers 18.40 18.17 
a3 19.07 18.70 Ton 7a 18.66 
g4 19.09 18.77 18.97 19.07 
C4 (19.85) 19.17 19.11 19.13+ 
c2 (19.99) 18.94 19.20 19.16 


* Local standards from SHP II. 
+ Variable ?. 
** Variable (mean of three measures). 


magnitudes derived from three photographic 
plates in each color. The faintest standards, how- 
ever, were only measurable on one of the three 
plates in each color. 


* Publication of the Goethe Link Observatory of Indiana 
University, No. 34. 
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Plate I. Photoelectrically measured standards are identified by lettered numbers and letters; most of them are i 
the vicinity of a, b and c. Numbers identify bright stars measured in field and listed in Table V. The photograph, taker 
with the 74-inch reflector at Radcliffe Observatory, reaches limiting magnitude in V. } 

3 
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te II. Bright stars measured in the center of the cluster are listed in Table II. South is at the top, west is at the 
left rather than on the right as has been the custom elsewhere in this series. 
Cc 


te I1I. Numbers identify faint stars measured in the NGC 458 cluster and adjoining field. Their magnitudes and 
; colors are listed in Tables III and IV, respectively. 


178 


3. Photographic colors and magnitudes. The 
central regions of NGC 458 are shown in Plate 
II, a short exposure taken in good seeing with 
the 74-inch Radcliffe reflector diaphragmed to 
44 inches. Most of the bright stars in the cluster, 
down to a limiting magnitude of about V = 18, 
could be measured on two such plates in each 
color. The derived magnitudes and colors are 
listed in Table II. In some cases where the stars 
were very central or crowded, eye estimates were 
made to obtain the magnitudes of these stars 
relative to the nearby stars with minimum error 
due to background light. 


TABLE II. BRIGHT STARS IN NGC 458—REGION I 
(Area = 2.501’) 


Star V B-V Star V B-V 
I 15.14 Sy 31 15.63 507" 
2 15.97 .10 32 15.97 .15 
3 16.49 85 33 16.32 -54 
4 15.82 — .04 34 16.90 .83* 
5 16.91 .82 35 17.50 — .09 
6 16.54 56 30\ 8 17453) 07, 
i 16.81 15 37 16.43 13 
8 16.61 55 38 17.82 —.09 
9 16.18 50 39 D752 == 202 
10 17.90 —.05 40 17.40 —.24* 
II 18.07 = 106" 41 17.58 02 
12 18.07 03 42 Ly 70 00 
13 15.85 85. 43 16.40 50* 
£ ap 15 = . 44 16.40 80* 
e 6 a eS 48 16 y S 
17 16.80 70* 47 16.64 83 
18 17.82 05 48 Wage = Pe 
19 18.07 OZ 49 16.81 95 
20 17.40 = EP 50 DAs Fe alle! 
21 17.72 (.O1) 51 15.70 17 
oa F7*32 — .31 52 8 32 54* 
eee Seren 
25 16.15 89* 

26 15.80 1.00* 

27 15.80 00* 

28 16.41 -58 A 
29t 15.81f - 487 “ 

30 16.60 .60* 


* Eye estimate. 
+ Variable. 


Plate III identifies the fainter cluster stars 
which were measured further to the edge of the 
cluster on longer exposure plates. All stars within 
the annulus were measured—one plate in each 
color. The magnitudes of these stars are listed 
in Table III. Plate III (Region III) also shows 
the fainter stars measured in the field; the colors 


and magnitudes for these stars, all of which stay 


inside the 1/5 radius, are given in Table IV. 
Finally, a sample of bright stars over the 
larger region IV shown in Plate I was made to 
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give representation of bright SMC members 
the color-magnitude diagram of the field. 
colors and magnitudes for stars included in th 
homogeneous bright sample are given in Table \ 

4. Color-magnitude diagram of NGC 458. F 
ure I shows the color-magnitude diagram for t 
cluster as plotted from the stars listed in Tab 
II and III. The stars measured in the innerm 
regions are designated by open circles. The st: 
measured further to the edge on longer expos 
plates are designated by filled circles. The ove 
lapping of the two symbols in the diagram give 
a check that the two samples j join sooty a wit 
no systematic méasuring differences. 

Because of the small areas sampled in th 
cluster, the diagram is practically free from fie 
stars infringing from the general field of the SM 
Above V = 18 mag., the stars in Figure I a 
almost exclusively representative of the clust 
Fainter than V = 18, there are only a few st 
(presumably including all those of intermedia 
color) which are interlopers from the field. 

There are three outstanding features of thi 
color-magnitude diagram: (1) the large num 
of giants; (2) the fact that these giants ha 
colors of only B—V=.5 to .g mag. rather th 


TABLE III. FAINT CLUSTER STARS—REGION II _ 
(Annulus Area = 0.901’) 


Star Vv Bye 


Star V B-V 
I 18.66 Ss} 31 18.29 
2 19.20 == GOR} 32 17.65 a 
3 18.71 a6 33 18.42 
4 18.07 = 00 34 19.05 = 
5 19.12 = .08 35 18.84 =a 
6 19.00 100 36 18.51 7 
7. (19.68)! °—).10 27 19.40 
8 18.84 Sei} 38 18.64 
9 19.18 + .06 39 18.30 
10 19.37 10) 40 18.87 
II 18.46 15 4I 18.96 = 
12 18.48 -00 42 18.05 = 
13 18.27 =.18 43 19.04 
14 18.42 =U 44 19.03 
15 19.69 B22) 45 19.18 
16 19.11 .08 46 19.57 
E/E STG -98 47 -19..53 
18 19.02 20) 48 18.69 ie 
19 19.1IT —.17 49 18.35 
20 18.76 -O1 50 18.92 
21 18.07 -04 51 19.04 = 
22 18.75 -83 52 19.51 
23 19.57 .00 53 18.86 =a 
24 18.49 -O1 54 19.78 a 
25 18.24 1.04 55 18.96 am 
26 16.27 .50 56 18.31 = 
27, 19.03 -O1 57 18.42 = 
28 19.14 On 58 18.39 I 
29 18.10 HK 


59 June 


ing redder; and (3) the discontinuity or break 
about I mag. between the bright end of the 
ain sequence and those evolved giants. Another 
ay of stating (1) and (2) is to say that the 
ertzsprung gap, or more specifically the region 
Cepheid instability, is unusually narrow and 
2||-defined in this cluster. 

Using m—M=t19.2 mag., as a temporary 
orking modulus for the SMC, yields luminosi- 
ss of M,=—2 or —3 for the brightest stars 
1 the NGC 458 main sequence. If the uncer- 
inties of age-dating a cluster like this by the 
shhonberg-Chandrasekhar main-sequence turn- 
f point are not prohibitive, we can conclude 
at this cluster, like the Pleiades, must be some- 
nere around 2X10’ years old (Sandage 1957; 
hanson and Mitchell 1958). 

5. Variable star in NGC 458. Star No. 29 in 
able I is variable. Because it appears so near 
e center of the cluster, an accurate light curve 
puld be difficult to derive. But the star has an 
nplitude at least greater than .5 mag., and from 
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the two plates obtained here on which it could 
be clearly seen, approximate mean magnitude 
and color were estimated. These estimates place 
the variable near the gap in the evolved giants. 
It is possible that others of that group are un- 
detected variables. 

However, this gap in the evolved giants in 
NGC 458 appears to be analogous to the Cepheid 
instability gap and the variable is probably a 
Cepheid. More accurate mean magnitude and 
color might well place the variable directly in 
the gap. 

If the variable is a Cepheid its period should 
be about 4 days, corresponding to its apparent 
magnitude in the SMC. This is a common period 
for SMC Cepheids and strengthens the previ- 
ously mentioned similarities between NGC 458 
and what is known of the LMC cluster, NGC 
1866, which contains half a dozen 3-day Cepheids. 

6. Color-magnitude diagram of the field. Figure 
2 shows the color-magnitude diagram of the field. 


NGC 458 


NO 2 Oe ee eet 


LOR 12s AGG 
Bay 


sure 1. Color-magnitude diagram of NGC 458. Open circles represent bright stars from center; filled circles represent 
fainter stars in outer cluster annulus. 


ei 
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NGC 458- Field 


20 


-4 -2. 0.2 4.6 8 10 12 14 (16 (82032 
B-V 
Figure 2. Color-magnitude diagram of field adjoining NGC 458. Circled points represent faint sample in region I 
uncircled points represent bright stars from sample over additional area of region IV, and sequence stars. 
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TABLE IV. HOMOGENEOUS SAMPLE IN FIELD—REGION III 
(Area = 7.200’) 

Star V B-V Star V B-V Star V B-V Star V B-V 
I 17.85 77. 59 ~=16.83 53 116 =: 118.59 -60 Pe Ud lG/ 56 
2 18.53 725 60 18.79 -60 117 19.18 aH 1Ke) 174 19.01 .62 
3 19.12 .08 61 18.40 -79 118 19.38 -65 175 18.38 — 09 
4 19.31 -41 62 19.99 = 2 y/ 119 19.07 -4I 176 18.81 .26 
5 16.82 1.36 63 17.88 = 05 120 18.90 S47] 177 18.25 —.22 
6- 19.32 EO 64 18.30 -09 121 16.67 +03 178 19.93 —.20 
7 19.25 -64 65 17.07 1.09 122 TOZ13) -75 179 18.77 —.25 
8 19.12 20 66 19.99 —-10 123 17.98 =.09 180 19.13 -59 
9 18.69 —.,10 67 19.74 —.17 124 17.74 1.16 181 19.65 — .06 
Io =. 18.88 = Gilt 68 18.82 al) 125 eto 72 oii 182 19.17 resis 
II 18.01 -9I 69 18.52 —.02 126 16.72 3.53 183 19.01 -56 
12 19.11 -56 70 19.86 «03 127 7/1 -59 184 19.14 -53 
13 ~20.51 730) 71 20.29 —.42 128 19.48 —.12 185 19.82 .16 
ta.) 1Sn79 83 72. 20372 = 02 120 neeelOn22 1.52 186 19.79 — .09 
15 16.99 1.28 73 19.98 —. 16 130 18.11 a2, 187 19.40 —.45 
16 18.73 .05 74 18.25 —.21 131 18.97 .62 188 20.41 — .80 
17 16.96 1.25 75 19.00 — .08 132 W757 —.13 189 19.18 —.26 
18 17.98 SSH 76 18.70 gil 133 18.92 -54 190 =. 20.04 —!,EA! 
19 16.51 1.24 77 18.19 =).09 134 19.07 -65 IQI 19.69 .03 

20 19). 74 -47 78 19.23 .61 135 18.76 -61 192 LGA ate) he 

21 19.17 —.50 79 19.04 = LO 136 18.22 — .06 193 18.24 .56 

22 18.29 — 25 80 17.18 -90 137 18.32 .88 194 19.64 —.21 
23 19.04 pS 81 19.46 = alii 138 18.81 -48 195 19.36 TO, 

24 19.46 82 82 17.28 caer) 139 19.20 -41 196 =: 19.49 rT 
25 19.89 12 83 19.57 — slits) 140 18.30 = 07, 197 18.85 .66 
26 19.30 61 84 19.16 .83 141 18.80 — .26 198 16.84 .80 

27 19.66 41 85 19.13 — .05 142 18.48 —)..06 199 19.58 -50 

28, 20.12 —.14 86 17.86 —.05 143 17.329 54 200 17.83 —.14 
29 19.11 —.12 87 19.81 rly 144 17.49 1.10 201 19.68 135 

Bon | 19n29 -66 88 19.63 .82 145 18.58 -O1 202 19.04 +55 

31 18.68 —.15 89 18.48 .38 146 19.92 .26 203 19.09 -50 
32 18.88 On 90 19.49 Ae) 147 19.51 — .30 204 18.36 = 3 

33 19.14 -97 gI 19.53 .00 148 17.71 -92 205 18.09 1.01 
34 19.90 rd 92 19.77 = (oy! 149 18.61 —.23 206 20.04 —.AT 
35 18.67 -00 93 19.40 51 150 19.29 = 5h 207 19.20 47 
36 18.19 ri O) 94 18.95 =.19 I5I 19.14 -66 208 18.09 —.II 
37 16.25 = Ode) 95 19.94 52 152 19.69 -48 209 17.42 =.12 
38 TOL. +04 96 18.88 ay fi 153 18.89 aves) 210 ~=16.80 72) 
49) 19.92 -06 97 + 19.90 —a20 E5450) 19:72 09 211 16.86 1.30 
40 19.23 -4I1 98 17.34 1.14 155 19.71 -47 212 18.23 —.14 
41 19.13 waz 09)= 19.20 wok 156 18.76 eeO 213 eto 7O -54 
2 19.29 -50 100 )=—-119..55 —.15 157 16.74 81 214 19.00 -51 
43 18.72 -64 IoI 18.34 -85 158 17.75 — .21 215 18.86 73 
44 19.58 -06 102 eee) =.19 159 17.46 120 ZION 10).28 -36 
45 18.41 84 103 7 0 —.19 160 17-55 1.16 217 14.57 1 West 1 
46 19.40 340 104 19.13 02 161 DSE27, Seat 218 18.32 Sell 
47 19.61 -57 105 18.48 -51 162 16.58 81 219 18.48 —.09 
48 19.49 On 106 19.53 22 163 17.03 1.00 220 19.00 -58 
Ay © 10).79 ——nOE TO 7a EO Ay) 08 164 =16..38 -82 221 18.00 -88 
Oe) 519273 Af} 108 19.93 Se y/ 165 18.44 —.15 222 19.30 —.52 
51 18.28 E03 109 17.98 = 5.09) 166 19.30 -52 223 18.80 —i<12 
BZ (18.25 -76 IOs tOs27. -58 167 19.29 -51 224 19.48 Syl 
53 19.10 .60 III 19.41 — .08 168 17.57 .48 225 LO. 77 —.27 
54 17.99 1.32 112 19.31 .50 169 19.52 ==! 21 226 15.66 — 04: 
55 18.37 63 113 19.76 24 170 19.08 as) 227 16.53 1.42 
56 §=17.99 =.09 TIA 919228 SHE) 171 19.20 -52 228) 19.03 -54 
57 16.72 9/3 II5 19.63 — .28 172 17.54 =—.15 229 17.20 1.26 
58 18.78 -78 


fhe circled points represent stars in Table [V— 
hat is, all stars within a circular area of 7.20 
quare minutes of arc. In order to define the 
wighter regions of the color-magnitude diagram 
or the field, the remaining uncircled points were 
neasured. They represent additional bright stars 
measured inside a larger annular area surround- 
ng the deep sample. The resulting color-magni- 


tude diagram is remarkably similar to the ones 
obtained for the SMC fields adjoining NGC 419 
and NGC 391. 

We may conclude that this type of color- 
magnitude diagram is accurately representative 
of the general background of the whole east wing 
of the SMC, and presumably characteristic of 
the entire SMC. 
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TABLE V. BRIGHT STARS IN FIELD—REGION IV 
(Annular Area = 18.010’) 

Star 4 B-V Star 4 B-V Star Vv B-V Star ¥ B-V 
I 5 iy fe .82 49) 45-09 -51 76 = 18.38 —.04 IZ. 17539 1.08 
2 16.61 1.27 40 17.58 -43 FATE, 17.20 -80 II4 i7-Z7, — Ff 
3 16.42 —.14 41 17.86 — .35 78 15-62 -66 115 17-58 —.i 
A E754: .87 42 18.10 —.19 719 17-41 -36 116 18.40 -. 
5 17327 —.16 43 17.20 —.,14 18.23 -61 117. 36512 = 
@>. 46:87 .87 44 18.74 — .36 81 17.09 ~I2 118.- 37287 
Fe a oe 45 16.42 —.18 82 "- TB S4i 4 1 eeg5) 119 16.87 I 
8 16.99 45 46 13.94 —.32 83 17-41 —.19 120 15-57, 

9 18.11 —.22 47 16.51 -g2 84 17.76 5 a 121 19.03 - 

10 = 18.20 —.II 48 18.61 —.22 85 17.86 I.15 122 18.65 _ 
II 18.42 —.16 49 17-99 —.42 86 16.54 Daa 123 37266 
12) 18-24 -93 50 18.51 —.16 87 =17.91 -95 124 18.21 _ 
13 17-57 —.17 51 17.91 — .30 88 18.20 Og 125 (88237 = 
14 16.72 —.19 52 18.24 1.05 89 15.87 —>~-09 126 18.32 - 
15 17.86 —_ 24 53 17.47 1.06 go 17-79 a9, 127. <i feas E 
16% “37-89 56 54 +#17.89 —.18 gI 18.35 -43 128 16.90 i” 
17 L742 229 55 18.02 —.19 g2 18.53 —.13 129 =6©18.12 — 
18 17.52 —.21 56 17-41 1.16 93 17.50 Sa 2 130 16.71 
19 18.04 —.15 57 16.69 32 94 16.97 =L19 131 18.51 

7s dames iy Ref 83 58 16.49 .69 95 17.87 —.19 132 18.66 

21 LAE 1.00 59 18.51 = 254 96 §=18.37 -34 133° 16272 = 

225 Taye 1.19 18,21 —.24 seamen ee (1.07) 134 17.50 

25 Ee gi .89 61 17.98 —.22 98 18.57 (.88) 135° igeqn 

24 + 16.49 65 62 G77 —.18 99 +=18.27 —.16 136 =: 18. 18 - 

25 18.07 (.91) 63 18.47 —.15., 100 16.94 -69 137 i842 = 

26 17.58 (1.25) 64 17.09 .87 101 16.93 1.29 138 17.50 - 

27 16.82 44 65 17.69 —=323 102 16.63 1.49 139 §=6.: 16.34 a 

28 17.49 1.00 66 17.85 -77 103 17.86 1.06 140 =15.95 a 

29 16.48 dble 67 18.01 -03 104 17.61 1.01 I4I 15-47 
30 =: 16.51 II 68 18.08 x27 105 17-77 —.18 142 18-25 = 

at 18.34 — .34 69. 37550 -70 106 =: 18.02 —.20 143° J i7-7e - 

32 18.19 —.II 70 16.97 1.24 107 17.99 Bi 144 16.67 is 

33 17.92 —.22 71 17.82 —.22 108 18.52 -76 145 15-65 = 

34 16.67 1.01 72 bis —08 109 17-35 —E36 146 15.38 

35 16.61 89 73 17.91 ee 10'S “GSAT — ee 147 17.69 
36 18.06 72 74 15.97 -05 Ill 16.90 +.82 148 17.42 

37 17.61 — 19 75 17.21 —.26 112 18.02 1.03 149 16.71 

38 16.90 1.19 


7. Conclusions. \t will be shown in the follow- 
ing papers of this series that NGC 330, another 
“blue” cluster in the SMC, confirms the appear- 
ance of the color-magnitude diagram of NGC 
458. Further analysis and comparisons will be 
made at that time. It is possible, however, to 
make two important conclusions, at this time. 

First, an earlier conclusion is strengthened, 
that the color-magnitude diagram of the general 
background of the SMC is made up of older 
stars, such as those present in the NGC 419 and 
361 clusters, plus a group of considerably younger 
stars. These younger stars comprising the field 
are now seen to be represented quite well gener- 
ally, and even in detail by the kind of stars just 
discussed in NGC 458. In fact, simply super- 
posing the two individual color-magnitude dia- 
grams of NGC 419 and NGC 458 reproduces the 
color-magnitude diagram of the field of the SMC 
very well. 

Second, because the color-magnitude diagram 
of NGC 458 is so different from clusters of the 
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same age in our own Galaxy, we conclude | 
recourse to the Vogt-Russell theorem (Arp 19 
that the chemical composition of these SMC s 
is different. Again, because we see examples - 
high metal content in young stars of our ov 
Galaxy, it is suggested that we are seeing 

other possible case in the SMC, that of 8 
stars with low metal content. Regardless of t 
details of the differences, however, it is no 
necessary to abandon the hitherto useful work 
hypothesis that all galaxies contain the sar 
kind of stars. /a 
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RECOVERY OF THE ANDROMEDIDS 


By GERALD S. HAWKINS, RICHARD B. SOUTHWORTH, ano FRANCIS STIENON 
Harvard College Observatory and Boston University 


Receiwed February 3, 1959 


Abstract. The spectacular displays of the Andromedid meteors in 1872 and 1885 have never been repeated and they have 
ten been listed as a “lost stream.’’ Precise photographic measurements with the Harvard Super-Schmidt cameras over 
|e period 1952-56 have disclosed 23 meteors that follow the orbit of Biela’s comet and represent the remnants of the 
arm. The individual orbits show evidence of perturbations caused by the earth and Jupiter. It is deduced that the orbit 
|| the Andromedid swarm now passes 0.03 a.u. above the earth’s orbit. The Andromedid meteors described here have 
sen perturbed in the plane of the comet orbit, a process similar to the dispersion observed in the sunward tail of the 


‘met Arend-Roland. 


Introduction. Meteor astronomy was founded 
a the observation of spectacular meteor dis- 
lays in the nineteenth century. The history of 
te Andromedid swarm is well known and has 
een adequately summarized (Olivier 1925; 
isher 1926; Lovell 1954). Biela’s comet, parent 
) the swarm, split into two parts that were last 
‘en in 1852. Prior to the demise, the comet suf- 
‘red perturbations by Jupiter which caused a 
togressive shift in the nodes of the orbit. Gradu- 
‘ly the comet orbit approached that of the 
urth; intersection occurred during the period 
800-50 and then the orbits separated again. At 
ast two meteor streams have been associated 
ith the comet. The first stream was observed in 
irly December (Dec. 6-10) at irregular intervals 
etween 1741 and 1847. We doubt whether the 
ite of occurrence of meteors ever exceeded 100 
er hour for this stream. Admittedly, Brandes 
‘corded a rate of 400 per hour in 1798, while 
‘avelling on top of a mail coach, but he applied 
dubious correction factor to allow for obscured 
ision (Benzenberg and Brandes 1800). The 
scond stream occurred in late November and 
as first observed on Nov. 29, 1850, by Heis. 
t was this second stream that produced intense 
leteor storms in the years 1872 and_1885, and 
aowers of diminishing intensity in 1892 and 
899. The observations thus indicated a periodic- 
y of 62 years. At each return, as summarized in 
able I, the shower occurred at an earlier calen- 
ar date, because of the regression of the node of 


the mean orbit. The date of maximum and the 
radiant position have been reduced to the stand- 
ard equinox of 1850. It was established beyond 
doubt that these meteors were moving in the 
orbit of Comet Biela, and, like the comet, were 
being subjected to perturbations by Jupiter. 

Observations of both streams in the present 
century have been very sparse, and by 1904 the 
rate had dropped to about 20 per hour. In 1940 
Dole reported an hourly rate of 30 on Nov. 15 
and Prentice a rate of five during the period 
Nov. 27 to Dec. 4. This anticlimax led many 
investigators to believe that the shower had suf- 
fered severe perturbations and has been thrown 
beyond the reach of the earth’s orbit. Thus by 
1954 we find the Andromedids listed as a “‘lost 
stream” (Lovell 1954). 

When meteor streams become dispersed it is 
essential to apply precise measuring techniques 
so that orbits can be computed for individual 
meteors. The early work of the Harvard 
Meteor Program (Whipple 1954) showed two 
meteors with similar orbits which were desig- 
nated as “Stream VI.” It became apparent that 
these meteors were possibly remnants of the 
Andromedid swarm and that a vestige of the 
stream was still in existence. This paper de- 
scribes a systematic search for the remnants of 
the great swarm from photographic records ob- 
tained with the Baker Super-Schmidt cameras. 

Observational methods and data. The Super- 
Schmidt cameras installed in New Mexico in 


TABLE I. THE ANDROMEDID SWARM 


Hourly Radiant 

Date rate R.A. Dec. Obs. Reference 
1850 Nov. 29 5? 15° +62° Heis Heis 1877 
1867 Nov. 30 7? ape +48° Zezioli Schiaparelli 1871 
1872 Nov. 27 3,000 25° +43° Many Fisher 1926 
1885 Nov. 27 15,000 25° +45° Many Fisher 1926 
1892 Nov. 23 6,000 Re +40° Many Fisher 1926 
1899 Nov. 24 150 he +42° Many Fisher 1926 
1904 Nov. 21 20 26° +44° Bohlin Bohlin 1905 
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1952 carried out a continuous photographic 
survey, interrupted only by periods of moonlight 
and adverse weather conditions. The records for 
the month of November have been systemati- 
cally searched each year from 1952 to 1956. Fig- 
ure Ia, which gives the total operational time for 
each calendar date, shows that the coverage is 
fairly uniform. All meteors with radiants located 
between RA 0° and 50°, and Dec. 0° and 50° were 
selected as possible Andromedids. These meteors 
are listed in Table II together with two meteors 
of Stream VI photographed with small cameras 
in 1950. The radiant positions are for the equi- 
nox of 1950. From the period of Comet Biela we 
would expect Andromedid meteors to have an 
observed velocity in the atmosphere of 20 km/ 
sec. A velocity histogram for the meteors, given 
in Figure 2, shows a pronounced grouping of 
velocities between I9 and 21 km/sec. For com- 
parison we studied the velocities of sporadic 
meteors that occurred at other months of the 
year, having the same elongations from the apex 
of the earth’s way as the meteors of November. 
The distribution curve obtained from these 
velocities, shown in Figure 2, may be taken as 
representative of sporadic meteors. We may 
therefore use velocity as a further criterion for 
identifying Andromedids; they are marked ‘‘A”’ 
in Table II. We note from the histogram, how- 
_ ever, that of these meteors selected on the basis 


FIGURE la 


FIGURE Ib 


Figure ta. Total number of observing: hours for each calendar date in No- 
vember for the years 1952 through 1956. 
1b. Mean ee rate of Andromedids for November 1952 through 


1956. 
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of velocity some 30 per cent may nevertheless f} 
sporadic and not’ related to the Andromedij 
swarm. : 


the month of November is shown in Figure 1 
the years 1952-56 inclusive. The meteors used ‘i 
this curve of activity are those marked ‘‘A”’ if 
Table II. We found activity in the stream 
year and no periodicity could be detected fr 
the observations. As indicated by Figure 
the rate of photographic meteors was approx} 
mately 0.17 per hour between November 2 
23, witha pronounced 1 increase on November 
From a comparison of photographic and vi 
observations we estimate that the maxim 
rate of 1.0 meteors per hour would correspon 
to a visual rate of 5. 7 

Meteors that were identified as possi 
Andromedids were reduced by one of th 
methods. A complete reduction with an ex 
correction for deceleration in the atmospher 
was carried out on the longer photograp 
trails, by the method of Whipple and Jacchi 
(1957). Other meteors were reduced by tk 
“short trail’? method (Hawkins and Southwo 
1957) or by the “‘graphical’’ method (McCro: 
1957). In reductions by the last two methods 
approximate correction was applied for atm 
pheric deceleration. iy 

Orbital elements for each meteor were com 


20 25 30 Nov 


59 June 


TABLE II, METEORS FALLING WITHIN THE 
LIMITS OF THE SEARCH 


Velocity 
Nov. Apparent radiant in atmos. 
Date UR R.A. ‘ Dec. km/sec 
1950 6.42 34.6° + 8.5° 20.2 A 
7.35 25.2 +27.7 20.8 A 
1952 7.09 EG +41.7 22.2 A 
Fotl 20.5 +28.2 20.9 A 
7.14 13.9 +32.2 20.3 A 
12.23 23.6 +31.7 20.3 A 
3 12.23 26.5 +33.0 10.8 
2) 12.26 24.5 +32.0 21.1A 
: 12.26 33-5 +10.0 24.2 
12.29 6.0 +12.0 15.1 
12.34 45-5 +22.8 4.7 
: 15.30 42.0 +30.5 22.0 A 
21.37 46.2 +18.0 19.3 A 
1953 217; 33.2 + 6.0 22.6A 
2.34 5-9 47.5 11.0 
3-29 43-4 +19.5 23.0 
7.46 37.2 +16.1 19.8 A 
13.25 29.9 ae BRS 16.1 
"1954 13.09 3.4 +10.6 16.1 
| 16.15 23.5 +31.0 14.5 
{ 17.12 42.4 +33.2 12.5 
¥ 17.18 18.6 +27.5 17.3 
| 18.11 37.8 +14.5 19.7A 
22.16 ek i203 12.2 
22.18 8.5 +44.5 38.7 
22.19 25.0 +39.2 19.6A 
23.32 36.5 +14.5 17.0 
30.31 36.1 +35.0 14.3 
1955 5.16 44.0 2a 38.3 
5.18 26.3 +21.0 61.3 
5.18 21.3 + 8.0 13.3 
10.21 37.8 +16.7 21.8A 
11.43 32.2 +29.6 19.6A 
11.21 45-7 + 5.3 17.0 
II.40 = 45.7 +29.5 19.5 A 
14.27 36.0 +10.2 18.4 
14.29 36.1 +26.0 19.0 A 
14.30 40.5 +35.0 21.6A 
1955 14.34 34.7 +38.2 20.2 A 
14.23 23.0 +47.6 17.6 
22.26 37.3 +34.0 18.2 
1956 5.27 32.8 +14.2 23.8 
6.24 4.8 +38.4 20.1 A 
6.26 22.3 +15.5 20.4 A 
7.36 34.8 +22.4 19.2 A 
9.31 30.3 + 3.9 13.0 
9.32 29.7 +20.8 20.8 A 


uted by the method of Whipple and Jacchia 
1957). The radiant position was corrected for 
iurnal aberration and zenith attraction and re- 
uctions were carried out for the equinox of the 
ear in which the meteor appeared. Table III 
ives the orbital elements for each meteor and 
idicates the reduction method employed. The 
robable errors of each element were deduced by 
erturbing the radiant and velocity and recom- 
uting the orbit. The method of reduction used 
ffects the size of the probable error, as indicated 
1 Table IV. The position of the node, 2, depends 
nly on the date of observation and has negligi- 
le error. The longitude of perihelion, z, is ac- 
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curately obtained by the first two methods, but 
the error in the semi-major axis, a, is large. The 
accuracy of all three methods is sufficient for 
statistical purposes such as the determination of 
mean elements. We have included in Table IV 
the error in Ty ,Tisserand’s criterion, computed 
for Jupiter. 


Percentage 


12 15 20 25 30 35 40 
Velocity km/sec. 


Figure 2. Velocity distribution of meteors observed during 
the months of November, 1952-1956. 


Discussion. The existence of an Andromedid 
stream in December is open to doubt. McCrosky 
and Posen (1959) have identified a minor stream 
called the 6 Arietids occurring between Dec. 8 
and 13. This stream has 2 = 259°2 which is the 
same as the node of Comet Biela in 1772. The 
other orbital elements, especially z and 7, do not 
agree with those of the comet, and the relation- 
ship between the meteors and the comet has not 
been definitely established. The 6 Arietids and 
other visual meteor streams observed during 
December are not, in the authors’ opinion, di- 
rectly related to the Andromedid swarm; they 
will not be discussed further in this paper. 

Good agreement exists between the mean 
orbital elements of the Andromedids and the 
last observed orbit of Comet Biela, as shown in 
Table V. It is instructive to study the varia- 
tions in the elements of Comet Biela at each 
return. We have plotted the five elements, a, e, 
1, O, and zw in Figure 3; in all elements progressive 
changes appear which are caused by perturba- 
tions from the planets, particularly Jupiter. 
Observations of the comet are available for 80 
years during which period the comet suffered no 
unusually large disturbance in its path. Since the 
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Nov. a 
Date Le Be a.u. e 
53 2ae7, 2.77 0.763 
56 6.24 5.02 0.829 
56 6.26 4.85 0.837 
50 6.42* 3.12 0.763 
52 7.09* 2.79 0.732 
52 Fula 2.59 0.708 
52 7.14f 3.09 0.738 
50 7-35" 3-59 0.789 
56 7.36 2.00 0.639 
53 7.46* 2.29 0.679 
56 9.32 3.56 0.789 
55 10.21 2.48 0.733 
55 11.40 1.75 0.647 
55 11.43 2.54 0.697 
52 12.23f 3.17 0.749 
52 12.26 3.40 0.781 
55 14.29 2.14 0.651 
55 14.30 2.10 0.678 
55 14.34 2.40 0.688 
52 15.30 2.43 0.723 
54 18.11 2.34 0.678 
52 PCN 2.52 0.695 
54 22.19 3.34 0.748 


* Complete reduction. 
+ Reduced by short trail method. 
Un-marked meteors reduced by graphical method. 


Andromedid swarm followed the orbit of Comet 
Biela very closely, we may presume that the 
swarm would suffer perturbations similar to 
those of the comet. It is therefore not unreason- 
_able to extrapolate the changes in the orbital 
elements of the comet from its last appearance 
in 1852 to the present day. The mean angular 
elements, 2, 2, and z, of the Andromedid meteors 
agree almost exactly with these extrapolated 
values (Table V). The agreement of the ascend- 
ing node is still better if we determine it from the 
time of maximum rate (Fig. 1b); then Q is 
230°5. If we presume that the present elements 
of the Andromedid swarm may be determined by 
such an extrapolation, we find that the orbit of 
the swarm now passes between the earth and the 
sun, as shown in Figure 4. D, the closest approach 
of the earth to the swarm, is given by the ex- 
pression 
D = RAQ sini (1) 


where AQ is the angular distance between the 
node of the swarm and the position of closest ap- 
proach. From the present values of the elements 
we find that the swarm is closest to the earth’s 
orbit when the sun’s longitude is 220°4; at that 
time D = 2.5 X 108 miles. 

It is interesting to draw the projection of the 
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with major axis in the plane of the comet orbit 
We presume the ends of the ellipse to lie below 


7 Q T g 

(deg) (deg) (deg) (a.u.) 
8 4 117.0 0.655 
on ae 88.6 0.860 
2.5 223.7 100.4 0.788 
3.3 43-4 108.7 0.741 
17.4 224.6 109.6 0.749 
8.4 224.6 108.9 0.758 
11.2 224.6 100.2 0.808 
8.7 224.4 106.8 0.756 
3.8 224.8 116.9 0.724 
0.0 224.7 113.5 0.735 
4.1 226.7 109.7 0.753 
0.3 226.9 124.3 0.663 
7.2 228.1 135.5 0.619 
73 ‘ 228.2 II1.0 0.769 
10.0 229.7 106.9 0.796 
10.4 229.8 107.3 0.789 
5-3 230s. 118.6 0.748 
Lee 231.1 128.9 0.674 
12.2 231.1 | 117.4 0.748 
8.4 232-2 128.9 0.674 
1.6 55-2 120.8 0.754 
0.1 59.0 121.9 0.767 
12.14 239.3 108.4 0.842 


plane of the swarm on the surface ‘‘abcd” 1 
Figure 4. This projection is the view that wo 
be seen by an observer on the sun, and is sho 
in Figure 5. The earth passes directly below 
swarm qn November 3 and cuts through 
plane of the orbit of the swarm on November 
Thus the cameras have recorded a maxim 
photographic rate at the time when the ea 
passes through the plane of the orbit. We ma 
infer that perturbations have acted on the swa 
to spread the members out in the original orbi 
plane. A dispersion such as this has been ob 
served many times when comets disintegrate, 
The latest notable example was Comet Are 
Roland 1957 which developed a narrow 
pointing towards the sun. Whipple (1957) int 
preted this sun-ward tail as illuminated deb: 
lying in the plane of the comet’s orbit. From tk 
curve of activity in Figure 1b we may infer the 
the earth first enters the Andromedid stream on 
November 2 and leaves it on November 22. Th 
observations are consistent with the assumptior 
that the cross section of the stream is an ellips 


the earth on November 23 and above the eart 
on October 14. 
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TABLE IV. PROBABLE ERRORS 
Velocity a 
Radiant (km/sec) (a.u.) La e i Ly 
+7’ ' 0.1 +0.1 +0°2 0.01 +0.1 +0.01 
+7’ +0.3 +0.3 +0°2 +0.02 +0.1 +0.04 
+3° +0.8 +0.8 + 4° +0.06 +1.6 +0.1 
TABLE V. MEAN ORBITAL ELEMENTS 
Reduction Elements 
method Year a e t . Ty 
Complete 
vidromedid Short trail 1950-55 2.90 0.732 7.5 (225.5) 108.2 0.539 
steors All 1950-56 2.88 0.728 6.3 228.1 113.5 0.564 
‘ment } Observational 1852 3153 0.756 12.6 247.3 110.5 0.486 
ela Extrapolation 1954 3.41 0.780 7.5 230.5 108.2 0.486 
: 


' From the cross-sectional diagram of the stream, 
‘igure 5, we can now see why the averaged values 
} semi-major axis a and eccentricity e of the 
‘ndromedid meteors do not agree with the 
trapolated values of Comet Biela. If the values 
| a and e were equal to the extrapolated values 
| Table V, the meteors would pass 2.5 X 108 
jiles above the earth on November 3. If the 
jeteors are to be detected at all by collision with 
+e earth, then they must be perturbed from the 
ean position of the swarm. If the perturbation 
'to move the particle in the original orbital 
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Figure 3. Orbital elements for Comet Biela and 
the Andromedids from 1772 to 1956. 


plane, then the elements a and e, and only these 
elements, will be changed. It is of interest to see 
whether such changes in e and a can be produced 
by the gravitational effect of Jupiter. A simple 
test is to compute Tisserand’s (1896) criterion, 
Tz, given by the expression 


2 COS 7 


yea va(I faa Ci 


i 


a (2) 


I 
a 
where P, is the period of Jupiter. This quantity 
should remain constant within the limits of ex- 
perimental error if Jupiter is the only perturbing 
influence. The straight lines in Figure 3 are 
drawn on the basis of a constant value of Ty. 
They describe accurately the secular perturba- 
tions of Comet Biela throughout its history from 
1772 to 1852. The extrapolation predicts, there- 
fore, a value of 7; = 0.486 for the present orbit 


Figure 4. Three-dimensional drawing of the orbits 
of Comet Biela and the earth. 
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Figure 5. Cross-sectional view of the Andeomedia stream taken in a 
plane perpendicular to the orbit of the earth. 


of the swarm. The mean value of 7; for present 
day Andromedids is 0.539 + 0.01 as given in 
Tables IV and V. This discrepancy shows that 
Jupiter alone is not responsible for the observed 
changes in e and a. 

How much of the original swarm now remains 
intact is uncertain. To take the most extreme 
point of view, we could assume that the original 
swarm has been dispersed uniformly throughout 
the area whose boundaries are shown in Figure 5. 
Thus if the original cross-sectional area of the 
swarm were s and the stream required ¢ days to 
cross the orbit of the earth, then we would ex- 
pect the number NV of meteors observed per hour 
to be reduced to a value NV’; that is, 


Nst 


N' = aarp (3) 


where s’ is the present cross section of the stream 
and ?’ is the number of days required for the 
stream to pass across the orbit of the earth. 
From the observed duration of the swarm in 
1875, we deduce that its thickness was 0.4 X 108 
miles and that its cross section s = 0.126 X 10” 
square miles. From the fact that thesswarm was 
not observed in 1878 we may deduce that ¢ is less 
than 100 days. From the cross section shown in 
Figure 5, s’ = 354 X 10” square miles, and 
t! = 6.6 X 365 days. From Eq. 3 we'deduce that 
N' = 0.2 visual meteors per hour. This value 
agrees quite well with: the average hourly rate 
found in the present series of measurements. 
It is therefore possible that the original swarm 
has been considerably dispersed during the last 
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Extrapolated Position 
of Andromedid Swarm 


one hundred ork by planetary perturbat 
and other causes. It is also possible, howe 
that perturbations produced a “‘coma’’ of mete 
from Comet Biela some time before the sw: 
appeared. ; 

We feel certain that the meteors described in! 
this paper represent the vestigial remains of the 
great Andromedid swarms of 1872 and 1885, 
that further research will permit a more deta 
account of the perturbations, both gravitati 
and non-gravitational, which have acted on 
meteors during the last one hundred years. — 

Acknowledgments. The research in this p 
was supported jointly by the Army, Navy 
Air Force under contract with the Massachu: 
Institute of Technology. 


REFERENCES 


Benzenberg, J. F. and Brandes, H. W. 1800, Sternschnu 
pen (Hamburg: Friedrich Perthes), 80. } 
Bohlin, K. 1905, Astromoniska Iakttagelser och Und 
sokningar (Stockholm), 8, No. 2. | 
Fisher, W. J. 1926, Proc. Nat. Acad. Sci. 12, 728. 
Hawkins, G. S. and Southworth, R. B. 1957, Son Mi 
Contr. ‘Ap. I, 207. 
Heis, E. 1877, Resultate (Miinster), pp. 31,159,175. 
Lovell, A. C. B. 1954, Meteor Astronomy (Oxford Tinive 
ity Press). 
McCrosky, R. E. 1957, Smithsonian Contr. Ap. 1, 215. 
McCrosky, R. E. and Posen, A. 1959, A. J. 64, 25. 4 
Olivier, C. P. 1925, Meteors (Baltimore: Williams a 
Wilkins). . 
Schiaparelli, J. V. 1871, Sternschnuppen (Stettin, Germ 
edition Be G. von Boguslavski), PP. 92, 100. f 
ae . 1896, Traté de Mécanique Celeste (Paris) 


Whine F. L. 1957, Sky and Telescope 16, 426. 

——. 1954, A. J. 59, 201. 

Whipple, F. L. and Jacchia, L. G. 1957, Smithsonian Com 
Ap. I, 183. 


59 June THE ASTRONOMICAL JOURNAL 189 
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| 

Abstract. The moments of the frequency distributions of proper motion, of radial velocity, and of apparent magnitude 
)| the G5 stars of the 6th magnitude have been combined analytically to yield the moments of the absolute magnitude 
(tribution. Comparison of these observed moments with those deduced from assumed combinations of Gaussian dis- 
butions indicates that 25 per cent of the G5 stars of this brightness are dwarfs. This result agrees well with the observed 


} 

Analyses of spectral-type distributions to ob- 
‘in the space density of stars require a knowl- 
+ ge of the relative numbers of giants and dwarfs 
+ a function of apparent magnitude. For the 
‘inter apparent magnitudes, m > 10, adequate 
ita usually can be obtained directly from the 
tectrum-luminosity classifications in the area of 
‘esky concerned. The number of stars of a given 
jectral type is generally sufficient to yield rea- 
inably sound statistical results. For the brighter 
jars, however, the data are insufficient unless a 
irge area of the sky is included. It becomes im- 
»rtant then to gather as much information re- 
rding the percentage of dwarfs as possible from 
‘study of stellar motions. This data, supple- 
-enting and supporting that derivable from an 
‘amination of the spectra of many stars, then 
\n be used effectively in studies of stellar space 
\stribution. 

Unfortunately, when the percentage of dwarf 
jars among those of a given spectral class as 
)termined by a study of stellar motions is com- 
‘red with that derived by direct luminosity 
jassification of the spectra, serious discrepancies 
‘e evident among the G5 stars. This has been 
‘dicated previously and is well known (Mc- 
‘uskey 1952). Some of the discrepancy is due to 
ie inclusion of all stars Go to G8 in the class 
\lled G. The percentage of dwarfs changes so 
pidly with spectral type over this range that 
‘ich a grouping leads to uncertain results for the 
urrower G5 class. The present investigation is 
/1 attempt to derive from a statistical study of 
‘ellar motions a more definitive value for the 
)rcentage of dwarf stars among those of spectral 
youp G5. Modern machine methods of tabula- 
jon and calculation make the analysis relatively 
sy. A byproduct of the present study is a group 
| programs for the IBM 650 digital computer 
‘hich will perform the moment calculations for 
ry desired data. These programs are available 
» other workers on request. The method used 
‘as been applied to the Gs stars of the 6th appar- 


jrcentage of dwarf stars of Class G5 obtained from spectral surveys in sampling fields in the Milky Way. 


ent visual magnitude with the results discussed 
in this paper. 

Method of analysis. The method of analysis 
rests upon an investigation by Kiewiet de Jonge 
(1948). The fundamental assumption is that the 
distribution function of the space velocities of the 
stars concerned is invariant over the volume of 
space in which the stars are embedded. Kiewiet 
de Jonge introduced a frequency function of 
stellar velocities which refers to an average dis- 
tribution of stars over the whole sky. Galactic 
concentration of the stars, or any other uneven 
characteristic of their distribution in the sky, is 
removed. 

Let p be the absolute radial velocity of a star, 
uncorrected for solar motion. Let g(p) denote the 
frequency function of p obtained by averaging the 
frequency distributions tabulated for equal 
standard areas. This function is obtained directly 
from the observational data. We assume that the 
g(p) distribution is normalized, that is 


J ewer =I 


Similarly let f(u) denote the distribution func- 
tion of the total proper motions, u. This function 
also is obtained by averaging the frequency 
functions tabulated for equal standard areas of 
the sky. This function is not normalized. 

The third frequency function needed is that 
for the apparent visual magnitude m. We denote 
this by 4(m) and require that it be normalized. 

Kiewiet de Jonge (1954) has shown that the 
first four moments of the frequency distribution 
in absolute magnitude, g(), are related to the 
moments of the distributions in yp, p, and m by 
the equations: 


M=m+5logu—5logp+ 6.870 (1) 
M? = —m?*+25 (log u)’—25(log p)?+3.88 (2) 


M?=m'+ 125 (log u)*—125 (log p)*—18.45 (3) 
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M‘= —m'+625 (log u)*— 625 (log p)4 
ae +581.7 (log p)?—6m? M? 
—150(M?+m") X[ (log p)?—0.155 ]+81.22 (4) 


The bars in these equations refer to moments 

about the means of the distributions involved. 
Suppose, now, that the luminosity function 

¢(M) for stars of a given narrow range in spectral 


class is 
Le (M — M,)? 
¢(M) = 1S exp | - CY | 
a = (Mg)? 
ee | =a | (5) 


where a is the fraction of dwarf stars among those 
of magnitude m, 1 — a is the fraction of giant 
stars, M,, o, and Ma, ca are the mean absolute 
magnitudes and dispersions for the giants and 
dwarfs, respectively. 

By use of the IBM 650 calculator, it is a rela- 
tively easy task to compute the moments about 
the mean for a theoretical g(M) for varying 
values of the parameter a. The method used here 
for determining the percentage of dwarf stars 
among the G5 stars of the 6th magnitude is to 
compare the moments M, M?, M* and M‘ derived 
by Equation (5) for various values of a with those 
obtained from the observational data by use of 
' (1) to (4). 

The observational data. For this study, the 
6th-magnitude G5 stars have been chosen. 
Radial velocities, proper motions and apparent 
visual magnitudes for the G5 stars with 5.51 
< m, < 6.50 were taken from the Bright Star 
Catalogue (Schlesinger and Jenkins 1940). Supple- 
mentary information concerning the probable 
errors in proper motion was taken from the 
General Catalogue of 33342 Stars (Boss 1937). 
Wherever the Bright Star Catalogue was deficient 
in radial velocity data or in parallax data, an 
attempt was made to fill the gap by using the 
General Catalogue of Stellar Radial Velocities 
(Wilson 1953) and the Catalogue of Stellar Paral- 
laxes (Jenkins 1952). The number of G5 stars 
included between the magnitude limits noted 
above is 417. All pertinent data for these was 
punched onto IBM cards. A program was devised 
for computing on the IBM 650 the total proper 
motion, the tangential velocity, and the space 


velocity for each star having the requisite data. | 


The resulting values were also tabulated on the 
data cards. No corrections for solar motion were 
made. 
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only 322 had measured radial velocities. 
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Table I exhibits the distribution of the 4 
stars in apparent magnitude. The last colur 
gives the normalized distribution ~(m)dm. 


TABLE I. DISTRIBUTION IN APPARENT 
VISUAL MAGNITUDE 

Relative No. 
my No v(m) dm] 
5-51-5.60 20 .04.796 
5.61-5.70 28 .06715 
5.71-5.80 28 .06715 
5.81-5.90 22 -05275 
5.91-6.00 46 -11030 
6.01-6.10 31 :07433 
6.11-6.20 49 -11750 
6.21-6.30 41 .09828 
6.31-6.40 81 .19428 
6.41-6.50 7p .17030 
Total 417 1.00000 


In order to obtain an averaged distribution ¢ 
proper motions for the whole sky, the latter wa 
divided into 94 equal areas according to the p 
published by Campbell (1928). The 417 
were then distributed over these areas. On 
assumption that antipodal areas have the sai 
weight, we have taken this weight to be 


W, = 417/47 (m1 a5 N2) (k = 1,2,°% 


where 7 is the number of stars at (a, 4), m2 is’ 
number of stars at (a+ 12, — 6). The 
n(k, ) is the total number of stars lying in 
kth pair of antipodal areas (a, 6) and (a + 
— 6) and having proper motions between u 
pw + du, we have 


f(u) 


for the uniform sky. (7 

Table II shows the adopted distribution 0 
total proper motions for 417 stars, corrected Di 
Smart’s (1928) method for a probable error ¢ 
+0"0046 per year. 

The observed radial velocities of the G5 st 
in the magnitude interval 5.51 to 6.50 were 
analyzed in much the same way as that deser 
for the proper motions. Weights were assigned 
the pairs of antipodal areas depending on f 
numbers of stars in each according to the formu 


47 
PS DS nk, pu) W; 
k=1 


Wi = 322/47 (m1 + m2) 


where m1 and m» are defined above. Of the 4 
stars for which proper motions were availabl 


Table III gives the adopted distribution fur 
tion of radial velocities. The observed frequen 
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TABLE II. ADOPTED PROPER MOTION 
DISTRIBUTION—EVEN SKY 
F(») wu F(x) wh f(u) BK F(u) 

jo oO 105 4.28 205 1.32 B05q0 e145 
mueoes 150. 3:88 210 1.30 310 2.75 
@me28.2 115 3.53 215 1.28 320 1.79 
Sigmpao-3) «120 °3,18 220 _1.25 330 1.48 
S@estea 125 -2.88 225 1.23 335 1.77 
45 31.5 130 2.58 230 1.20 340 3.70 
° 30.5 135 2.43 235 1.15 360 1.92 
Gee25-3 140 2.28 240 I.10 395 2.01 
HO) 17-5 145 2.13 245. 1.07 400 0.81 
(5 14.7 150 1.98 250 1.05 405 2.00 
HO 12.6 155 1.83 255 1.03 430 0.81 
n5 11.3 160 1.69 260 1.00 465 1.11 
(0 10.1 165 1.59 265 0.97 485 177, 
N5 9-3. %170 1.50 270 0.95 500 2.00 
‘o 8.45 175 1.48 275 0.90 15 stars with uy >.5 
'5 7-70 180 1.45 280 0.85 
0 689 185 1.42 285 0.78 
'5 6.34 190 1.40 290 0.75 
10, 5-79 195 |1.38 295 0.77 
5 5-24 200 1.35 300 0.88 
i 4.68 


inction was smoothed graphically and the 
noothed normalized values are shown in the 
ible The integration under the smoothed g(p) 
irve, required for the normalization, was carried 
1t by Simpson’s rule on the IBM 650. No cor- 
‘ction in the distribution was made for accidental 
‘ror in the radial velocities. 

Figures 1 and 2 show graphically the observed 
|stributions in total proper motion and in radial 
»eed. The vertical line through each point in 
igure I indicates the natural uncertainty in the 
umber represented by the point. If N is the 
umber in a given interval, N? is the natural 
facertainty. Each vertical line is 2? in length. 
urve I in Figure 1 indicates the original 
noothed distribution in yw. Curve II is the dis- 
‘bution corrected for an accidental error in pu 
f +0%0046 per year. Curves III and IV are 


aried curves used to test the sensitivity of the 


TABLE III. ADOPTED DISTRIBUTION OF RADIAL 


: 
| VELOCITIES—NORMALIZED 


p 
(km/sec) AO) (km/sec) (o) 

fo) -O41901 65 -0008717 
5 -039660 70 -0006849 
10 -034991 75 .0006226 
15 .029138 80 -0005603 
20 .022476 85 .0005603, 
25 .015877 90 -0004981 
30 -009775 95 .0003736 
35 -006600 100 -0003113 
40 -005043 105 .0003113 
45 -003985 110 .0003113 
50 -002864 115 .0002490 
55 -001930 120 -0002490 
.001245 125 -0001 868 
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moment calculations to reasonable changes in 
the f(u) distribution. 

Analysis of the data. Kiewiet de Jonge (1954) 
has shown that the sampling errors arising even 
in samples of 1000 make the third and fourth 
moments of the absolute magnitude distribution 
computed by Equations (3) and (4) very un- 
certain. In the case of the Ao stars of magnitudes 
6 to 7, which he investigated, the probable errors 
in the third and fourth logarithmic moments of 
the radial-velocity distribution amounted to 20 
per cent and 32 per cent, respectively, of the 
moments themselves. For the third and fourth 
logarithmic moments of the proper-motion dis- 
tribution, the corresponding figures were 15 per 
cent and 18 per cent. On the basis of a sample of 
1000 stars, the third and fourth moments of the 
absolute-magnitude distribution were shown to 
have probable errors in excess of +3 mag. The 


probable error in M? for the same size sample was 
found to be 40.3 mag., and that in M, +.06 mag. 
In view of the large uncertainties to be expected 
in the third and fourth moments, our principal 
conclusions will be drawn from the first and 
second moments of the absolute-magnitude 
distribution. 

The logarithmic moments for the smoothed 
total proper-motion and radial-velocity distri- 
butions for the G5 stars were computed on the 
IBM 650. Simpson’s rule for numerical integra- 
tion was used. For comparison, a direct-moment 
calculation using the data tabulated in classes 
was also made. The results are shown in Table 
IV. Entries in the table are in every case mo- 
ments about the mean. For the f(y) distribution, 
the adopted values of Table II were used. 

By means of (1) to (4), the moments of the 
absolute-magnitude distribution of the G5 stars 
were found to be those listed in the last line of 
Table IV. 

An examination of Figure 1 indicates that there 
is considerable leeway in drawing the smooth 
curve through the points. The effect on the mo- 
ments in log » that would be caused by drawing 
other curves than that adopted (Curve I1) is 
exhibited in Table V. These varied curves are 
numbered J, III, IV in Figure 1. Table V shows 
that the logarithmic moments are not very 
sensitive to small but possible variations in the 
drawing of the f(u) curve. 

By direct calculation with the IBM 650 the 
moments about the mean have been obtained for 
a hypothetical absolute-magnitude distribution 
as defined by (5). The parameters used are from 


30 


Frequency 
$ 


10 


Proper Motion ¢.001"/yr.) } 
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Figure 1. Frequency distribution of total proper motions. Roman numerals designate distributions as follows: F 
(dashed curve) Original smoothed distribution; II. Distribution corrected for accidental error; III and IV. Arbitraril? 
varied curves to test sensitivity of results. Lines through the points indicate natural uncertainties of the data. 


’ Keenan and Morgan (Hynek 1951) and Nassau 
and MacRae (1949): 
-M, =+0.2 
9 =t08 


Mz =+ Bp! 
oe = 10.7 


for the giants and dwarfs respectively. Table VI 
gives the first four moments about the mean for 
¢(M) as defined by (5) for various values of a, 
the fraction of dwarf stars involved. 

Direct interpolation of the observéd value for 
M from Table IV between appropriate values of 
ain Table VI yields a = 0.23. Similar interpola- 
tion for M? yields a = 0.28. These two values 
are in good agreement and indicate that about 25 
per cent of the G5 stars with 5.51 < m, < 6.50 
are dwarfs. 

As pointed out above, the higher-order mo- 
ments are practically useless in such an analysis. 
In the present case the values of M/* and M* taken 
from Table IV lie considerably outside the range 


of M? and M‘ of Table VI. Since M? is definitely | 


positive for the observed data, about all one can 
say is that a strong preponderance of giant stars 
is indicated. 


.and M? from this observed y(M) distributio 


It should be noted in passing that the values ; 
a derived from M? is not unique. The observer 


value of M? fits into Table VI just as well a 
a = 0.73 as at a = 0.28. However, the high e 
value is inconsistent with that derived from J 
and hence may be rejected. Zz 

These results for a may be compared with sim 
ilar values deduced from a distribution in ab: 
lute magnitude of those G5 stars for which M 
be computed directly. Of the 6th magnitude 
stars, 203 have absolute magnitudes given” 
the Bright Star Catalogue, or computable fre 
more recent data. For these, a distribution 4 
cording to M in o. 5 magnitude intervals w 
made. This is shown in Table VII. The momer 
about the mean for this distribution compar 
with the moments deduced from a study of # 
motions are shown in Table VIII. The third at 
fourth moments are very discordant and ine 
cate, in general, the fluctuations which are causet 
by sampling. If we interpolate the values of 4 


between the entries of Table VI, we find ae 
0.29 and 0.17 respectively, the average beim 
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TABLE IV. MOMENTS ABOUT THE 


Method Distribution log u 
Direct Calc. F(u) — 1.3066 
| Simpson’s Rule F(u) —1,3020 
: 
| log p 
| Direct Cale. g(p) 1.0467 
Simpson’s Rule g(p) 1.0336 
m 
Direct Calc. h(m) 6.1152 
M 
Formulae 1-4 Absolute Mag. 1.3069 


= 0.23. This is not greatly different from the 
ilue found by the analysis of the motions. 
These results are in good accord with the per- 
tage of dwarf stars found for all stars of class 
5 in the Bright Star Catalogue for which abso- 
‘te magnitudes are available. If we arbitrarily 
tablish the dividing line between giants and 
warfs at M, =+ 2.4, we find that 89 of the 
+6 stars have M, > 2.4. This is 24 per cent of 
ie whole. 


| TABLE V. VARIATION IN MOMENTS ABOUT THE 
MEAN—f(y) DISTRIBUTION 


| Curve loge (log x)? (log )3 (log 1)é 

| I. 1.2711 0.29008 0.08752 0.32569 

etl —1.2942 0.30378 0.09981 0.29566 
Ill —1.2675 0.30877 0.05287 0.34747 
IV —1.2635 0.29611 0.04919 0.33198 


TABLE VI. MOMENTS ABOUT THE MEAN— 
THEORETICAL ¢(M) DISTRIBUTION 


| a mM M2 Ms mM 
0 +0.2 0.64 o 1.229 
2 +1.180 4.452 10.942 61.076 
4 +2.160 6.342 5.118 58.781 
6 +3.140 6.312 —6.176 59.716 
8 +4.120 4.362 — 11.647 62.845 
1.0 +5.1 0.49 {o) 0.721 
| TABLE VII. DISTRIBUTION OF OBSERVED M 
FOR 6TH MAG. G5 STARS 
M No. M No. 
<= 3 (0) +2.01+2.50 4 
—2.99 — 2:50 4 +2.51 + 3.00 I0 
—2.49 — 2.00 3 +3.01 + 3.50 12 
—1.99 — 1.50 Oo +3.51 + 4.00 6 
1.40% 91,00; 3 +4.01 +450 6 
—0.99 —0.50 5 +4.51+5.00 6 
—0.49 — 0.00 8 +5.01 +5.50 9 
+0.01 +0.50 26 +5.51 +600 3 
+0.51 + 1.00 62 +6.01+ 6.50 I 
+1.01 + 1.50 18 +6.51 +7.00 5 
+1.51 + 2.00 12 27.0 ° 
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MEAN—6TH MAG. G5 STARS 
og x)? og 1) ‘(log x)4 
0.29756 0.09191 0.30896 
0.32020 0.08625 0.38254 
(log p)? (log p)* (log p)* 
0.22603 —0.17900 0.16600 
6 0.26300 — 0.16698 0.42438 
m ms m4 
0.0776 —0.0110 0.0124 
Mw Me Mt 
5.2324 13.1927 119.586 
TABLE VIII. COMPARISON OF MOMENTS 
FROM TWO DISTRIBUTIONS 
M mM Mi Me 
From Motions 1.307 5.232 13.193 119.586 
From MM Distr. 1.604 3.926 5.247 50.386 


TABLE IX. G5 STARS IN LF REGIONS 


(133 sq. deg.) 
Luminosity Class 


III IV Vv 
Mp No. % No. % No. % 
7.25 267. (o) 2 (0) Tiere a3 
7.75 (o) Ce) 014.0 2 100 
8.25 2 100 O10 fe) co) 
8.75 Aan O7, o oO 2 33 
9.25 9) 42 6 29 6 29 
9.75 16 57 Suats 7 25 
10.25 28 865 IS oh Fai 
10.75 35) 40 14 16 38 44 
11.25 ye 2) 16 12 59 45 
11.75 96 8646 Lig Bs oes 103 49 
12.25 149° 57 EG 99 38 


TABLE X. PERCENTAGE OF DWARF G5 STARS 
AS FUNCTION OF MAGNITUDE 


Mp NG Mp % 
6 25 10 45 
7 28 II 54 
8 33 12 65 
9 «38 


In the foregoing analysis no account has been 
taken of the subgiant stars (Class IV) of spectral 
group G5. Trials of a tri-modal distribution, 
¢(M), to include these stars did not yield results 
as consistent as those derived on the basis of the 
bi-modal theoretical distribution. The intermedi- 
ate Class IV stars in this analysis have been in- 
cluded among the dwarfs since their mean abso- 
lute magnitude is +3.4 according to Keenan and 
Morgan (1951). 

The fainter G5 stars. The spectral data accumu- 
lated for studies of the stellar luminosity function 
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Frequency 


pitkm/sec) 


Figure 2. Frequency distribution of radial velocities. 


at the Warner and Swasey Observatory (Mc- 
Cuskey and Seyfert 1947) indicate the trend in 
the percentage of dwarf G5 stars with increasing 
apparent magnitude. Nine of these LF regions 
cover an area of 133 square degrees and are 
spaced along the galactic equator from / = 12° 
to 1 = 182°. Table IX gives the numbers of stars 
and the percentages in luminosity classes III, IV 
and V for the G5 stars in these regions. If the 
percentages for'classes IV and V, the dwarf stars, 
are combined so that the results are comparable 
to those already derived for the brighter stars, 
we find the trend shown in Table X. The per- 
centage of dwarf stars increases by a factor of 2.6 
from the 6th to the 12th magnitude. The data 
presented in Tables [X and X refer to apparent 
photographic magnitudes. 

The results presented in Table X differ some- 
what from those deduced by other investigators. 
Oort (1936), for example, finds 42 per cent dwarf 
stars among those of classes G4—G6 for latitudes 
less than 30° and for m, = 11.5 to 13.4. This 
value contrasts with the 71 per.cent dwarf stars 
of class G5 found by van de Kamp and Vyssotsky 
(1937) for m, = 10.2 and B=0° to +10°. 
Vyssotsky and Williams (1948) found 45 per cent 
dwarfs among this class of stars for m, = 10.5 
to 11.4 and for all galactic latitudes combined. 
It is apparent that considerable scatter exists in 
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these determinations. This is evident also in 
percentages exhibited in Table IX. 

It is a pleasure to record our indebtedness 
the staff of the Case Computing Center for 
advice and facilities made available to us. 
IBM 650 programs are recorded in the reports Of 
the Computing Center and are available upon 
request. The Case Research Fund has generou: 
supported the computational work involved. 
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NOTICE 


DEFINITION OF THE NEW I.A.U. SYSTEM OF GALACTIC COORDINATES 


| 

At the Dublin General Assembly of the International Astronomical Union in 1955, Sub-Commission 
b of the Union was appointed “‘to investigate the desirability of a revision of the galactic pole and 
the zero of galactic longitude”. At the Moscow General Assembly in 1958 the conclusions and 
ommendations of this Sub-Commission were reported and discussed in a joint meeting of Com- 
ssions 33 (Structure and Dynamics of the Galaxy), 33b, and 40 (Radio Astronomy). The General 
‘sembly subsequently passed a resolution proposed by Commissions 33 and 40 which, among other 
‘ms, contained the following recommendations: 


“. . . (c) That Commission 33b be authorized to define the exact values of the coordinates of 
| pole and of the zero of longitude immediately after the final reduction of the relevant observations 
finished. 
(d) That Commission 33b be charged with the communication of these values to the members of 
e 1.A.U. and to all other interested institutions and individuals’. 
Commission 33b has now completed the final reduction of the observations referred to in item (c), 


id it has defined the following values of the relevant quantities: 


@ The new north galactic pole lies in the direction 
| 


Q =12549™ 6 = +27°4 (equinox 1950.0). 


‘2 position angle 


| (2) The new zero of longitude is the great semicircle originating at the new north galactic pole at 
. 6 = 123° with respect to the equatorial pole for 1950.0. 
\ 


] (3) Longitude increases from 0° to 360°. The sense is such that, on the galactic equator, increasing 
lactic longitude corresponds to increasing right ascension. Latitude increases from —90° through 
‘to +090° at the new north galactic pole. The system is therefore similar to the Ohlsson system 
coordinates. 

The above quantities are to be regarded as exact so that the new galactic coordinates may be 
puted to any desired accuracy in terms of right ascension and declination for the equinox 1950.0. 


| Other useful values (restricted accuracy). The quantities which follow have been computed from 
ie definition but are given only to an accuracy of the nearest integer in the last digit. 


@ The 1900.0 values corresponding to those in the definition 
| New Pole 
a-= 1254676, 6 = +27°40’ (1900.0) 
: New Longitude Zero at position angle 6 = 123°04’ (1900.0) 
©) The old galactic coordinates of the new pole 
P= 347°7, bt = +88°51 
(c) Position of the point of zero longitude and latitude (new system), I’ = 0, b = 0. 


In equatorial coordinates 


= 17"30%3, 6 = —28°54’ (1900.0) 
=a 42 44 6 = —28°55’ (1950.0) 


In old galactic coordinates 


J = 327°609, bt = —1°40 
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for galactic longitude and latitude respectively. Commission 33b eee chat, during the tran 

period, the symbols /!, b! should be used for the old system and /", 6! for the new one. The 
mission also strongly recommends that, apart from these superscripts, it should be made quite « 
whether the galactic coordinates used in any publication are based on the old or the new “yet 


longitude na through the galactic center. The pole position was based primarily on ee ial 
tion of neutral hydrogen in the inner parts of the Galaxy (R < 7 kpc, assuming Rsun = BS < 
the most precise evidence on the galactic center came from the location of the radio source S. 
tarius A, which is presumed to be the galactic nucleus. It remains a problem for future investigat 
to determine to what extent more exact positions of the layer of neutral hydrogen and the gals 
nucleus deviate from the quantities adopted for the present revision. However, it may be expe 
that these deviations are small and will not give rise to further revisions of the galactic coordi: 
system in the near future. 4 

The new method of specification of the zero of longitude in terms of its position angle has | ie 
chosen because this angle is one of the quantities used directly in conversion formulae from equat Ke 
to galactic coordinates. ; 

The equatorial plane of the new coordinate system must of necessity pass through the sun. 
a fortunate circumstance that, within the observational uncertainty, both the sun and Sagittarit iv 
lie in the mean plane of the Galaxy as determined from the hydrogen observations. If the sun | 
not been so placed, points in the mean plane would not lie on the galactic equator. 

A report will soon be published giving details of the Commission’s investigations. 

Conversion tables. In accordance with a decision taken at the Moscow Assembly, conversion + is 
from equatorial into galactic coordinates and vice versa, and from the old and new galactic coordit 
systems into each other, will soon be published under the supervision of Commission 33b. a 


The members of Sub-Commission 33 


A. BLAAUW 
C. S..GUMo ae 
- J. L. PAWSEY | 
25th March 1959 j G. WESTERHOU 


